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Abstract 

Discussed is a model of collective and internal degrees of freedom with 
kinematics based on affine group and its subgroups. The main novelty in 
comparison with the previous attempts of this kind is that it is not only 
kinematics but also dynamics that is afnnely-invariant. The relationship 
with the dynamics of integrable one-dimensional lattices is discussed. It 
is shown that afnnely-invariant geodetic models may encode the dynamics 
of something like elastic vibrations. 

Keywords: collective modes, affine invariance, integrable lattices, nonlinear 
elasticity. 

Introduction 

In some of our earlier papers including rather old ones [TTJ [501 [2D HH H31 GUI 
(MIMSlllSllHlllilSllMllMll^lSnilSllMllMll^ we have discussed 
the concept of affmely- rigid body, i.e., continuous, discrete, or simply finite 
system of material points subject to such constraints that all affine relations 
between its elements are frozen during any admissible motion. For example, all 
material straight lines remain straight lines in the course of evolution, and their 
parallelism is also a constant, non-violated property. Unlike this, the metrical 
features, like distances and angles, need not be preserved. In other words, such 
a body is restricted in its behaviour to rigid translations, rigid rotations, and 
homogeneous deformations. Models of this kind may be successfully applied in 
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a very wide spectrum of physical problems like nuclear dynamics [5] (droplet 
model of the atomic nuclei), molecular vibrations, macroscopic elasticity [T5lfT6l 
Sli3llllia[51[51[51[5Sl[Sni[Slll5](in situations when the length of excited 
waves is comparable with the size of the body) , in the theory of microstructured 
bodies [53J (micromorphic continua), in geophysics [H [TT] (the theory of the 
shape of Earth), and even in large-scale astrophysics (vibrating stars, vibrating 
concentrations of the cosmic substratum, like galaxies or concentrations of the 
interstellar dust). 

From the purely mathematical point of view such a model provides an in- 
teresting example of a system with the group-theoretical background of the 
geometry of degrees of freedom [5J [UJ ESI HS1 EH HH1 ES] ■ It is an affine general- 
ization of the usual rigid top with the orthogonal group replaced by the linear 
one (isometries replaced by affine transformations). Let us remind also that 
there is an interesting formulation of the general non-constrained continuum 
mechanics based on the infinite-dimensional " Lie group" of all diffcomorphisms 
or volume-preserving diffeomorphisms [TJ ^ [3J [T3J HU [33J 131] ■ This theory 
is rather complicated (although heuristically very fruitful) because of serious 
mathematical problems with infinite-dimensional groups. The mechanics of an 
affinely-rigid body is a simple compromise between rigid-body mechanics and 
such a continuum theory, because admitting deformative degrees of freedom it 
is simultaneously based on the finite-dimensional framework. 

Let us stress, however, that, in spite of its non-questionable physical applica- 
bility and formally interesting features, the referred mechanics of affinely-rigid 
body is in a sense disappointing from the point of view of the mathematical 
theory of Lie group motivated systems. The point is that in the latter the- 
ory it is not only kinematics (finally, geometry of degrees of freedom) but also 
dynamics that is ruled by the underlying group. Due to the isotropy of the 
physical space, Lagrangian of a free rigid top, i.e., its kinetic energy, is invariant 
under all left regular translations (all spatial rotations); the same is valid, of 
course, for the resulting equations of motion (Euler equations). If the material 
structure of the top is isotropic (spherical inertial tensor), then the model is 
also invariant under right regular translations. When formulating the theory of 
ideal incompressible fluids in terms of the group of all volume-preserving diffeo- 
morphisms, one obtains an infinite-dimensional Hamiltonian system invariant 
under all right regular translations. This is due to the fact that in the usual 
Euler description of fluid its Lagrangian coordinates are not very essential, and 
the fluid particles have a rather limited individuality. Summarizing, in these 
theories one deals with Lagrangians or Hamiltonians based on left-, right-, or 
even two-side invariant metric tensors on the Lie group used as a configuration 
space. It is never the case in the above-quoted model of affinely-rigid body. 
This brings about the question as to the hypothetic affine counterpart of left- 
and right-invariant geodetic models on the orthogonal group and their potential 
perturbations. This interest is at least academically motivated. But at the same 
time, from the physical point of view, such models look rather esoteric. In any 
case, the previously mentioned applications of affine collective modes are dy- 
namically well-established, because they are based on the d'Alembert principle 
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in theory of constrained systems. There are, however, some indications that 
physical applicability is not a priori excluded. 

This problem has to do with the very philosophy of the origin of collective 
and internal degrees of freedom. We say that a "large" system of material points 
(continuous, denumerable, or just finite admitted) has collective modes when 
there exists a "small" number of parameters q 1 , . . . , q n that are dynamically rel- 
evant, i.e., satisfy an approximately autonomous system of evolution equations, 
if for our purposes the kincmatical information about the system, encoded in 
them, is sufficient, and (very important!) if they depend on individual particles 
in a non-local way. The latter means that positions and velocities of all particles 
enter the (f-variables on essentially equal footing, with the same strength, order 
of magnitude, so to speak. This is, of course, a rough, qualitative introduction 
of the term, but there is no place here to develop a rigorous mathematical de- 
scription. As a mathematical model we can realize some quotient manifolds of 
multiparticle state spaces or their submanifolds (e.g., representatives of cosets). 
On the contrary, internal degrees of freedom are described in terms of fibre bun- 
dles over the physical space, space-time, or the configuration space. They give 
an account of phenomena which are either essentially non-extended in space, 
or perhaps cannot be described in terms of composed systems because their 
spatial details are unapproachable to our experimental abilities. For example, 
from the point of view of contemporary science, spin systems seem to be based 
on essentially internal quantities [301 171) . In any case, spin media do not look 
like the Cosserat continuous limit of discrete systems of molecular "gyroscopes" . 
The latter model works successfully in the theory of Van der Waals crystals and 
granular media. 

Apparently, the most natural and intuitive origin of collective modes, e.g., of 
some microstructure variables, is based on the mechanism of constraints and the 
d'Alembert principle. Collective motion is then "large", whereas non-collective 
one is "small" and merely reduced to some vibrations about the appropriate 
constraint submanifold. The collective kinetic energy, i.e., dynamical metric 
element, is obtained from the restriction of the total one to the constraints 
surface (the first fundamental quadratic form). This corresponds to the classical 
relationship between kinetic energy and inertia [2J [71 8, 69J. In this case, as a 
rule, the collective kinetic energy is invariant under a proper subgroup of a group 
underlying geometry of the constraints submanifold. But one can also realize 
another mechanism, namely, such one that the hidden non-collective motion is 
just large, and that the emerging collective modes have to do with the averaged 
behaviour of hidden modes, i.e., with the time dependence of some relatively 
slowly- varying mean values. Then it is quite natural to expect that the collective 
Lagrangian will be based on a kinetic energy whose underlying dynamical metric 
tensor will be non-interpretable in terms of the restriction of the usual multi- 
particle metric tensor of the kinetic energy to the constraints manifold (i.e., 
to the first fundamental form of constraints). Similarly, equations of motion 
need not be derivable from the usual d'Alembert principle based on the original 
spatial metric. Therefore, the relationship between kinetic energy and inertia 
may become rather non-classical, to some extent exotic in comparison with the 
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usual requirements (cf., e.g., the discussion by Capriz and Trimarco [TJIHIEI]). In 
such situations the only reasonable procedure is to postulate the kinetic term of 
the Lagrangian on the basis of some natural and physically justified postulates. 
Let us mention two examples from the two completely opposite scales of the 
physical phenomena, namely, the atomic nuclei and vibrating-rotating stars (by 
the way, the neutron stars are in a sense exotic and gigantic nuclei with Z = 
and enormous A). As objects more close to the Earth one can think, e.g., 
kinetic bodies as discussed by Capriz, and various non-standard microstructure 
elements like gas bubbles, voids, and defects in solids [71 [51 US]- Though bubbles 
and voids can be hardly treated as constrained pieces of a substance or systems 
of material points. 

Situation is even much more complicated, when one deals with essentially 
internal degrees of freedom, like, e.g., spin systems [71] . Then, although we 
have some guiding hints from the theory of extended systems, any choice of 
Lagrangian, Hamiltonian, or equations of motion is based on some rather hy- 
pothetic postulates, first of all, on certain invariance requirements. 

There is also another point worth of mentioning. Namely, usually in varia- 
tional theories of analytical mechanics, Lagrangian consists of the kinetic and 
potential parts. The first one has to do with inertia, constraints, metric struc- 
ture, whereas the other one describes true interactions. But even in traditional 
problems of analytical mechanics there are approaches where the structure of 
interactions is encoded in an appropriate metric structure, i.e., in a kind of 
kinetic term. There is a well-known example, namely, the Jacobi-Maupertuis 
variational principle. If ds is the usual metric (arc) element of the configuration 
space, and V is the potential energy, then one uses a modified metric [2] 



where E denotes a fixed energy value. This is so-called isoenergetic dynamics, 
based on the homogeneous "Lagrangian" 



A denoting an arbitrary parameter (not time). This variational principle, based 
on the metric element day, gives trajectories with the energy value E, but 
without the time-dependence. There are also spatiotemporal forms of this prin- 
ciple, where the time variable occurs as one of coordinates q l , and there is no 
restriction to the fixed energy value. 

In a slightly different context, in certain problems we will follow the idea of 
encoding the interaction structure in an appropriately postulated kinetic energy 
form, i.e., metric tensor on the configuration space. 

As mentioned, we concentrate below on models with kinematics (and dy- 
namics) ruled by the linear group GL(ro,R), or, more rigorously, affine group 
GAf(n, R) (physically n — 2,3). Of course, the usual rigid body in n dimensions 
is ruled by SO(n,R), or, if translations are taken into account, by the isome- 
try group E(n,R) = SO(n,R) x s R™. But there are also other possibilities of 



day = V E — V ds, 
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finite-dimensional collective modes, e.g., SL(n,R) (or SL(n,R) x s I n ), i.e., in- 
compressible affinely-rigid body, or, just conversely, the Weyl group R + SO(n, R) 
generated by rotations and translations (the shape of the body is preserved, but 
not necessarily its size). In some future we are going to investigate systems ruled 
by the projective group in n dimensions, Pr(n,R) ~ SL(n + 1,R), cf., e.g., [55] , 
This is quite a natural extension of affinely-rigid body, when the system of ma- 
terial straight-lines is preserved but their parallelism may be violated. Another 
interesting model would be given by the Euclidean-conformal group CO(n,R). 
Let us mention that there was also some very interesting attempt by unjustly 
forgotten German physicist Westpfahl [72], who invented the idea of using the 
unitary group U(3) as a basis for collective modes in three dimensions, quite 
independently of later applications of unitary symmetry in elementary particle 
physics. 

Finally, it is quite often so that the complexification idea leads to physically 
interesting results. It is not excluded that complexifying the physical space R n 
to C™ and replacing the real groups GL(n, R), U(n) by GL(n, C) we could obtain 
some interesting models of collective or internal degrees of freedom [7511751175] . 
The idea is particularly tempting, because GL(n, R) and U(n) are two different 
(and qualitatively opposite) real forms of the same complex group GL(n,C). 
But, of course, such exotic ideas are rather far from realization and they are 
mentioned here only because of their obvious conceptual relationship with the 
usual and generalized models of affinely-rigid bodies. 

The group space is a particular model of systems with kinematics and dy- 
namics ruled by a Lie group. In general, the microstructure or collective con- 
figuration space (the manifold M. in the sense of Capriz book [7j and related 
papers) is a homogeneous quotient space G/H. Here G is a fundamental group 
of the model, and H is an appropriate subgroup of G, not necessarily normal 
one, thus G/H need not inherit the group structure from G [7J[8]. 

1 Dynamical systems based on Lie groups 

Dynamical systems based on Lie groups and their homogeneous spaces are 
widely used as a model of internal and collective degrees of freedom [7J [8] [9] 
[TOl [3T1 [32] . They present also interest by themselves from the purely mathe- 
matical point of view. They are realistic and quite often they possess rigorous 
analytical solutions in terms of special functions and power series; this is prob- 
ably due to the analytical structure of Lie groups. The first step of analysis is 
the theory of left- and right-invariant geodetic systems, when the Lagrangian 
and total energy are identical with the kinetic energy expression based on an 
appropriate Riemannian structure of G. 

For simplicity let us use the language of linear groups; by the way, nonlinear 
groups are exceptional in applications, and the most known examples are the 
universal covering groups GL(n,R), SL(n, R) of the indicated linear groups. 
For any curve R 3 t \— > g(t) S G its tangent vectors g{t) S T g t t \G may be 
transported to the Lie algebra G' = T e G with the help of right or left g(f) -1 - 
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translations, resulting in quantities f2(f) := g^g^)^ 1 , Ct(t) := gft)^ 1 ^^). In 
this way the tangent and cotangent bundles TG, T*G may be, in two canonical 
ways, identified with the Cartesian products: TG ~ G x G", T*G ~ G x G'*. 
It is clear that the left and right regular translations g <— > Lk(g) = kg, g i— > 
Rk{g) = 3& transform quasi- velocities either according to the adjoint rule or 
trivially: 

L fe : I— > Adfeil = fcr^fc !li->4 

i? fe : fi^ft, f2 i— > Adfc-if2 = k^Clk. 

Left-invariant geodetic systems on G are based on kinetic energies, which 
are quadratic forms of ft with constant coefficients. If G is non-Abelian, then fl 
is a non-holonomic quasi-vclocity and the corresponding Riemannian structure 
on G is curved. Similarly, right-invariant kinetic energies are quadratic forms 
of Q with constant coefficients. As a canonical example of left-invariant sys- 
tems we can realize the free rigid body in n dimensions, G = SO(n,M) (if we 
neglect translational motion). If the rigid body is spherical (its inertial tensor 
is completely degenerate), then T is also right-invariant, and the underlying 
metric tensor on G is proportional to the Killing tensor. Such a pattern may be 
followed in all semisimple Lie groups [21 [33 O [42] . Quite a different example 
is provided by the theory of the ideal fluids |2j . The configuration space is iden- 
tified with SDiff R 3 — the infinite-dimensional group of all volume-preserving 
diffeomorphisms of R 3 (provided that we discuss the physical three-dimensional 
case). If we admitted the fluid to be compressible, we would have to use the 
full group Diff R 3 of all diffeomorphisms. The functional of kinetic energy is in- 
variant under right regular translations in SDiff R 3 . What concerns left regular 
translations, it is invariant only under the six-dimensional isometry group of R 3 . 
The reason for this relatively poor left-hand-side invariance is that the kinetic 
energy expression depends in an essential way on the spatial metric tensor. At 
the same time, from the point of view of the material space, the particles of fluid 
have a rather limited individuality, and that is why the kinetic energy form of 
incompressible fluid is invariant under the huge group of sufficiently smooth and 
volume-preserving "permutations" of particles, i.e., under SDiff R 3 . 

In some of our earlier papers [T^l [5111 1211 1221 IH3 ISHl ESI E31 EH SHI SSI 
E3[54l[5g[60l[Ml[63l[64l[S6]we discussed the object called " affinely-rigid body", 
i.e., the system of material points "rigid" in the sense of affine geometry, i.e., all 
afflne relationship between constituents being kept fixed during any admissible 
motion. Such a model is geometrically interesting in itself and has a wide range 
of applications in macroscopic elasticity, mechanics of micromorphic continua 
with internal degrees of freedom, molecular vibrations, nuclear dynamics, vi- 
brations of astrophysical objects, and the theory of the shape of Earth [JJ ITT] . 
Analytically, the configuration space of n-dimensional affinely-rigid body may 
be identified with the semi-direct product GL(rt,R) x s R", or simply GL(rt,R) 
when we neglect translational degrees of freedom. 

The kinetic energy of an extended affinely-rigid body in Euclidean space 
may be calculated in the usual way, by summation of kinetic energies of its 
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constituents. Velocity vectors are squared with the use of the fixed metric 
tensor of the physical space. The resulting metric tensor of the configuration 
space is flat, and it is not invariant either under left or right regular translations, 
except two subgroups isomorphic with the Euclidean group SO(n,R) x s R n . 
Because of this the resulting geodetic model, although kinematically based on 
the group manifold, dynamically is incompatible with it. Besides, it is physically 
non-realistic and useless, because geodetics are straight lines in L(n,R) x s R ra , 
therefore, in certain directions the body would suffer a non-limited extension or 
squeezing. It is impossible to avoid such non-physical catastrophic phenomena 
without introducing some potential term. 

The very taste and mathematical machinery of systems with group-manifold 
degrees of freedom consist in the invariance of geodetic models under the to- 
tal group of regular translations. This motivates the search for left- or right- 
invariant kinetic energies, i.e., Riemannian structures on GL(n,R) x s R™ or 
GL(n,R). The first step is purely mathematical: the very construction and 
some primary analysis. Later on some hypotheses are formulated concerning 
the physical applicability of such apparently exotic "non-d'Alembertian" mod- 



2 Kinematics and Poisson brackets 

Let us remind briefly the basic ideas concerning the extended affinely-rigid body 
in a flat Euclidean space [ill [13 HSl SHI E21 E31 EH ESI E01 - It is convenient 
to use the standard terms of continuum mechanics, although the model applies 
also to discrete or finite systems of material points (provided there exist at 
least n + 1 material points in rt-dimcnsional space) . Two Euclidean spaces are 
used, namely, the material space (N, U, 77) and the physical space (M, V, g); the 
symbols N, M denote the underlying sets, U and V are their linear spaces 
of translations, and -q e U* <£> U* , g € V* <X> V* are metric tensors. We put 
dim A = dimM = n. The points of N are labels of material points. The 
configuration space Q of affinely-rigid body in M is given by AQ.(N,M), i.e., 
the manifold of affine isomorphisms of A onto M. Obviously, it is an open 
submanifold of AfI(A, M) — the affine space of all affine mappings of TV into 
M (including non-invertible ones). In some configuration $ 6 Q the material 
point a £ N occupies the spatial position x = <5(a) E M. The co-moving, i.e., 
Lagrangian, mass distribution within the body will be described by the constant 
(time-independent) positive measure /i on N; it may be <5-like (concentrated 
at single points), continuous with respect to the Lebesgue measure, or mixed. 
Cartesian (Lagrange) coordinates a K in N are chosen in such a way that their 
origin is placed at the centre of mass C, i.e., 



The manifold Afl(JV,M) may be identified with the Cartesian product Mx 
L1(U,V), where Ll(U, V) denotes the manifold of all linear isomorphism of U 



els. 
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onto V; it is an open submanifold of the linear space L({7, V). The first factor 
refers to translational motion, i.e., to the centre of mass position x = The 
linear part of <&, ip = L[<f>] = 6 IA(U,V), describes the relative (internal) 
motion. Analytically, when some Cartesian coordinates in M are used, motion 
is described by the dependence of Euler (current) coordinates on Lagrangian 
(material) ones and on the time variable: 

$(t,o)* =<p i K (t)a K +x i (t). 

In practical calculations it is often technically convenient, although may be 
geometrically misleading, to identify both U and V with 1" and Q with semi- 
direct product GAf(n,R) ~ GL(ra,R) x s R™. Another natural model of Q is Mx 
F(V), where F(V) denotes the manifold of all linear frames in V. By the way, 
F(V) as a model of internal (relative-motion) degrees of freedom is essentially 
identical with LI(C7, V) if we put U = R™ and use the natural isomorphism 
between linear mappings tp € LI(R™, V) and co-moving frames e G F(V) frozen 
into the body and attached at the centre of mass. This must be done when 
the body is infinitesimal and the relative motion is replaced by the dynamics 
of essentially internal degrees of freedom. Then R™ becomes the micromaterial 
space of internal motion. 

Inertia of the body is described by two constant quantities, namely, the total 
mass and the second-order moment of internal inertia J € U <£> U , i.e., 

dn(a) 7 J KL := I a K a L dn{a) 

N JN 

(cf., e.g., 0HligilHlllSlliaiS3IS2EaEniE21); J is symmetric and positively- 
definite. 

Summing up the kinetic energies of constituents, 



1 f 3$* d& 



one obtains: 



rr rr , rr m dx l dx l 1 dip 1 A dip 1 B t ab ri x 

the symbols "tr" and "int" refer, obviously, to the translational and internal 
(relative) terms. 

The phase space of our system may be identified with the manifold P := 
M x L1(U,V) x V* x L(V,U). The factor V* refers to translational canonical 
momentum, whereas L(V, U) to the internal one, in the sense of the obvious 
pairing between tt e L(V,U) and generalized internal velocity £ £ L(U,V): 
(n, £) = Tr(7r • £) = Tr(^ • tt). Cartesian coordinates in M generate parametriza- 
tion pi, p A i of canonical momenta. For Lagrangians of the form L = T—V(x, ip) 
Legendre transformation 



dx j A dip 3 b T RA 

Vi=mg i:j — , p ^ i = g ij ——J (2) 
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leads to the following kinetic term of the Hamiltonian: 

T = ^ gt3ptP] + l9 tl P A t P B 1 JAB, 

where, obviously, g lJ are components of the reciprocal contravariant metric of 
g, and J e U* <£> U* is reciprocal to J, Jac'J CB — $A B , do not confuse it with J 
with the ^-lowered indices. This kinetic term (and its underlying flat metric on 
Q) is invariant under Abelian additive translations in Q — Mx Ll(U, V); those 
in the second term are meant in the sense 



Ll(U, V) 3 if ^ <p + a, aeL(U,V). (3) 

Therefore, without the interaction term (for L = T), the Hamiltonian genera- 
tors pi, p i are constants of motion. However, as mentioned above, such geode- 
tic models for deformable bodies are physically non-interesting, because they 
predict unlimited expansion, contraction, and passing through singular config- 
urations with det ip = 0. The latter, although non-acceptable in continuum 
mechanics, may be to some extent admissible in mechanics of discrete bodies. If 
we once decide that the internal configuration space is given by IA(U, V), then 
the above transformation group is only local. At the same time, even for purely 
geodetic systems, as mentioned, there is no invariance under geometrically in- 
teresting affinc groups of left or right affine regular translations in Q. Even if, 
at the present stage, models with affinely-invariant kinetic energy might seem 
rather academic, they present some interest at least from the purely mathe- 
matical point of view. Besides, some physical applications seem to be possible 
in hydrodynamics, astrophysics, nuclear dynamics, and in certain elastic prob- 
lems. It is interesting that even without any genuine interactions, on the purely 
geodetic level such models may predict bounded and stable elastic vibrations 
of incompressible bodies. It is so as if the interaction was encoded in the very 
kinetic energy, i.e., configuration metric, so as it is, e.g., in Jacobi-Maupertuis 
variational principle. To formulate such models we must introduce and partially 
remind certain geometric objects. 

Affine velocity in laboratory representation, i.e., expressed in terms of space- 
fixed frames, is defined as 

0:=^L<n (,-)*,. 

The corresponding co-moving object, related to the body-fixed frame, is given 
by 

6- ,„-i^ cT mi 6-4 (, n -i\ a d( P l B 

n . = ip -euu), n B = (<p ) 4 — 

Obviously, ft — ip&ip -1 , = ip" 1 a& A b (v? 1 ) B j- These are Lie-algebraic ob- 
jects corresponding to the structure of Q as the group space of a Lie group. 
They provide an affine counterpart of the rigid-body angular velocities, and in 
fact reduce to them when ip is confined to the manifold of isometries of (U, rj) 
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onto (V, g); then they become skew-symmetric respectively with respect to 77 or 
9- 

The object f2 may be represented in terms of continua as a gradient of the 
Euler velocity field, namely, the material point passing the hxed spatial point y 
has the translational velocity: 

E v{yy = — + rti{yi-xl), 

i.e., simply £l l jy J in the instantaneous rest frame of the centre of mass, placed 
also at the instantaneous position of this centre in M. Similarly, dp> % s/dt has to 
do with the gradient of the Lagrange velocity field, because the instantaneous 
velocity of the a-th particle (a e N) is given by 

v ' dt dt 

(concerning the standard concepts of continuum mechanics consult, e.g., [15[ 
UniUS])- In certain problems it is also convenient to express the centre of mass 
translational velocity v % — dx l /dt in co-moving terms, i.e., v A = (<^ _1 ) A iV % ■ 

It is very convenient to introduce the canonical affine spin, also in two 
representations, the spatial and co-moving ones £ G L(V), £ £ L(U). In 
terms of coordinates they are given by the following formulas: £*j = p> % aV j> 
t A B = p A i p> i B - As previously, S = p±<p~\ S%- = ^ A t A B (p- 1 ) B r 

They are purely Hamiltonian quantities defined on the phase space; with- 
out any precisely defined Lagrangian or Hamiltonian we cannot relate them to 
generalized velocities. It is seen, however, that they are dual objects to affine 
velocities, i.e., they are non-holonomic canonical momenta conjugate to them 
in the sense of following pairing: 

(S, fi) = (S, n) := Tr(EO) = Tr(±Q) = p A i v i A , 

where v % a are generalized velocities of internal (relative) motion. This canonical 
isomorphism between Lie algebras GL(U)' = L([/), GL(V)' = L(V) and their 
duals simplifies remarkably all formulas and considerations. 

It is clear that quantities T?j are Hamiltonian generators of GL(V) acting 
on 1A(U, V) through the left translations: 

p ^ Asp, tpe Ll(U, V), A e GL(V). (4) 

Similarly, T, A b generate right regular translations in the internal configuration 
space: 

tp^ipB, p£Ll(U,V), B G GL(C/). (5) 

In continuum mechanics these mappings are referred to, respectively, as 
spatial and material transformations; in this case they include rotations and 
homogeneous deformations. Obviously, to use correctly such terms we must 
be given metric tensors in V and U. Then the g-antisymmetric part of E and 
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the ^-antisymmetric part of E generate, respectively, spatial and material rigid 
rotations; the symmetric parts generate deformations. 

The doubled antisymmetric parts are referred to as spin S and vorticity V 

[12, 

S*,' = E*j - g ik g jm Z m k , V A B = t A B - v AC VBD t D c . (6) 

Attention! There is an easy mistake possibility: if motion is not metrically- 
rigid, then V is not a co-moving representation of S, i.e., 

S*j * ^aV a b {p- 1 ) B r 

Just as translational velocity, the canonical linear momentum may be ex- 
pressed in co- moving terms according to the following rule: pa = Pity 1 a- 

The objects f2 and E are invariant under material transformations, but the 
spatial action of A E GL(V) transforms them according to the adjoint rule, 
i.e., £1 i ► AflA^ 1 , E i— ► AYiA^ 1 . On the contrary, Cl and E are invariant 
under GL(V) but experience the inverse adjoint rule under B E GL(U), i.e., 
fl i— > B~ 1 ClB, E i— > B~ 1 Y,B. This formally agrees with formulas for systems 
with configuration spaces identical with Lie groups, but one must stress that 
there are some subtle differences due to the fact that LI(J7, V) is not a Lie 
group (may be identified with it, but there is an infinity mutually equivalent 
identifications). 

The translational or orbital affine momentum with respect to some point 
O G M is defined as follows: 

where x % are Cartesian coordinates of the O-radius vector of the current position 
of the centre of mass in M. The total affine momentum with respect to O is 
given by 

i{oy 3 ^Acoft + sy 

A(O) and I(O) depend explicitly on the choice of O. Unlike this, E is objective 
(in a fixed Galilean reference frame). There is a complete analogy with the 
properties of angular momentum, the doubled ^-antisymmetric part of the above 
objects. The quantity I(O) is a Hamiltonian generator of the group of affine 
transformations of M preserving O (O-centred affine subgroup). 

Poisson brackets of E-quantities follow directly from the standard ones for 
x l , pi, tp l A : P A i- The non- vanishing ones are simply given by the structure 
constants of linear group, 

{E^EM = 5\Y, k 3 - S k jY; i i, {^ 3 ,t A B } = 0, 

{E^B, E C £i} = S C B^ A D — 5 A ]j'E C b 

(similarly for A, /). There are also non- vanishing Poisson brackets related to 
the left or right affine groups GAf(M), GAf(iV). Here belong the above ones 
and besides, those involving linear momenta, 

{t A B ,p c } = 5 A cf>B, {I l 3 ,Pk} = {V 3 ,Pk\ = 5 l kP r 
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If F is any function depending only on the configurations variables, then, obvi- 
ously 

BF BF BF 

{F,* i} = ^ ^ {F, A*j} = {F,t* B} = 

Geometric meaning of the last formulas is clear, because the differential opera- 
tors used on their right-hand sides are identical with vector fields on Q gener- 
ating the action of one-parameter subgroups of GAf(M) and GAf(iV). 

As mentioned, the above Poisson brackets follow directly from the standard 
definition [lEEUM] 

{ ' 1 ; ~ dq a d Pa d Pa dq a ' 

where q a are generalized coordinates and p a are their conjugate canonical mo- 
menta. In our model q a are given by x l , <p l a, and p a by pi, p A i. In applications 
it is sufficient to remember that {q a ,q^} = 0, { Pa ,Pp} — 0, {q a ,p@} — S a p, 
that Poisson bracket is bilinear (over constant reals K), skew-symmetric, i.e., 
{F,G}= -{G,F}, satisfies the Jacobi identity {{F,G},H} + {{G,H},F} + 
{{H, F}, G} = 0, and finally that 

k 

{F, H(Gi, . . . , Gfc)} = ^2 H, p (Gi, . . . , Gk){F, G p }, 
P =i 

where commas before indices denote the partial derivatives. The formerly- 
quoted Poisson brackets together with the above rules are sufficient for all cal- 
culations concerning equations of motion and their analysis. 

To define a non-dissipative (Hamiltonian) dynamical model, we must be 
given some Lagrangian L{q,q), perform the Legendre transformation, p a = 
dL/dq a , invert it, i.e., solve with respect to generalized velocities q a , and sub- 
stitute the result to the energy function E = q a dL/dq a — L. In this way one 
obtains the Hamilton function H(q,p). Equations of motion may be then ex- 
pressed in terms of Poisson brackets, 

£-<'■->• 

where F runs over some finite family of basic functions, e.g., (pi,Y, l j,x l ,(p l A), 
(Pa, ^ A b, x i , ip % a), or something else. The basic dynamical laws are given by the 
balance equations for the linear momentum and affine spin either in laboratory 
or co-moving representation (one could use equivalently the linear momentum 
and the total affine momentum, however, the previous choice is more conve- 
nient). The procedure based on Poisson brackets and canonical formalism is 
very often more easy and computationally less embarrassing than the one di- 
rectly using the Euler-Lagrange equations. 

Remark: Legendre transformation may be also expressed in terms of non- 
holonomic objects, moreover, this is often more convenient and effective than 
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the use of generalized velocities. Expressing Lagrangian in terms of (V, fi*j) or 
(v A , £i A b) instead of (p l a), we can describe the Legendre transformation as 
follows: 

dL . dL dL ~, dL 



dv v J W ™ 0**' 

When dealing with the Hamiltonian form of equations of motion, we need 
often Poisson brackets involving deformation tensors and certain by-products of 
the inertial tensor, like, e.g., the Eulerian quadrupole of the mass distribution. 

Obviously, for systems with affine degrees of freedom the Green and Cauchy 
deformation tensors G G U* <£> U*, C G V* <g> V* are respectively given by 
the following expressions: G = (p*g, C — (p^ 1 ) Tj, i.e., in analytical terms 

Gab = gij(fi l A<P 3 B, CV, = 77^5 (<P -1 ) A i ( ( y 5_1 ) B j- Their inverses G € U ®U, 
C G V (g> V are defined by G ac G C b — S A b, C ik C k j = 5 i j , and one must be 
careful to avoid mistaking G AB , G 4 -? with -q AC rj BD G C d-, g lk g jl Cki- Therefore, 
the usual convention of the upper- and lower-case indices may be misleading. 
Analytically, G AB = (ip' 1 ) A t (ip- 1 ) B 3 g zj , C ij = ^ A p j BV AB - When there is 
no deformation, i.e., p 6 LI([7, 77; V,g), then G = 77, C — g. The corresponding 
deformation measures vanishing in the non-deformed state, i.e., Lagrange and 
Cauchy deformation tensors E £ U* <£> £/* , e G <8> V * are given by (see, e.g., 

[mug): 

£:=i(G-77), e:=i(s-C). 

One uses also their contravariant versions .E" 45 , e* 5 '; unlike G AB , G u they are 
defined via the 77- and g-raising of indices. 

Remark: G is independent of 77 and may be defined even if the material space 
is purely affine, amorphous. Similarly, G is independent of g and is well-defined 
even if the physical space is metric-free. Therefore, the literally meant term 
"deformation" is better expressed by E, e than G, G. However, in many formu- 
las G, G are more natural and convenient. Deformation tensors behave under 
the action of isometries in a very peculiar way, namely, for any A G 0(V,g), 
B G 0(U, 77), we have: 

G[A(p)kl = G[<p]kl, GYpB]kl — G[ip]cD B C K B D L , 

C[Ap]^ = C[<pU (A- 1 ) \ (A- 1 ) \, C\pB\ l3 = C[dy. 

By the way, the last two formulas are valid for any A £ GL(V), B G GL(C7). 
The first two equations (invariance rules) imply the Poisson-bracket rules 

{G KL ,S i j } = 0, {Cij,V A B } = 0, 

and similarly for Ekl, ey- 

Deformation invariants arc important mechanical quantities. They are scalar 
measures of deformation, basic stretchings, which do not contain any informa- 
tion concerning the orientation of deformation (its principal axes) in the physical 
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or material space. They may be chosen in various ways, but in an n-dimensional 
space exactly n of them may be functionally independent. The particular choice 
of n basic invariants depends on the considered problem and on the computa- 
tional details. When non-specified, they will be denoted by K a , a = l,n. Let 
us define mixed tensors G £ U <g> U* , C £ V <g> V* 7 E £ U ® U* , e € V <E> V* , 
namely, 



G A 



B ■' 



V AC G 



CB, 



C l 



9 ik C kj , 



E- 



AC 



E, 



CB, 



9 ik e kj - 



A class of possible and geometrically natural choices of K, a is given by the 
following expressions: Tr(G fe ), Tr(C fc ), Tr(£ fc ), Tr(e fe ), k = l,n. In certain 
problems it is convenient to use the following eigenequations: 



det G b 



X5 A , 



0, det C 



0. 



det 



\5 J 



0, det [e*j - \S l j] 



0. 



These are n-th order algebraic (polynomial) equations with respect to A. Their 
solutions provide one of possible choices of basic invariants. Another, very con- 
venient one is given by coefficients at X p , p = 0, (n — 1) [HI EI] (the coefficient 
at A" is standard and equals one). Deformation invariants are non-sensitive 
with respect to spatial and material isometries, i.e., for any A £ 0(V,g), B £ 
0(U,rj) we have JC a [A(pB] = JC a [(fi]- This implies the obvious Poisson brackets: 
{lC a ,S* 3 } = {lC a ,V A B } = 0. 

In certain computational problems, but also in theoretical analysis, it is 
very convenient to use quantities Q a = y/X^, where A a are solutions of the 
above eigenequations, or q a — lnQ a (i.e., Q a = exp(g a )). The eigenvalues of C 
equal (Xa)- 1 = {Q a y 2 = exp(~2q a ). 

Any function F on the configuration space which depends on (p only through 
the deformation invariants is doubly isotropic, i.e., satisfies F(AipB) = F(ip) for 
any A £ 0(V, g), B £ 0(U, rj), <p S LI(t/, rj; V, g). All such functions have vanish- 
ing Poisson brackets with spin and vorticity, i.e., {F, S l j} — {F, V A b} = 0. In 
certain formulas we need the spatial inertial quadrupole, J[f] ab = f a Kf b ' lJ kl ■ 
It is related to J KL just as C is to -q. When the body is inertially isotropic, 
J[tp] becomes proportional to the inverse Cauchy deformation tensor. Unlike the 
co-moving internal tensor J £ U ' ®U, J [if] £ V ® V is configuration-dependent, 
thus variable in time. 



3 Traditional d'Alembert model 

At least for the comparison with more exotic (although geometrically and per- 
haps physically interesting) suggestions we must start with a brief reporting and 
extension of the traditional model based on the d Alembert principle. As shown 
in [13 gi SS EH ED , Lagrangians of the form L = T—V(x,ip) with T given 
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by ((T|) lead to the following dynamical laws: 



d Pi _ dV dE*j _ t , dV _ , 4 

~ ~ ^ ~dT- Umlj i-VAg^-Urn^ j + Qj, (7j 

expressed in terms of Cartesian coordinate systems. This is the balance for 
fundamental Hamiltonian generators. It becomes a closed dynamical system 
when considered together with the Legendre transformation |(5J) or its equivalent 
description 

dx j 

Pi = mgij—j-, E*,- = g jk tt m J[<p] m% . (8) 

Substituting these expressions to the dynamical balance for pi, S l j one obtains 
some reformulation of the Euler-Lagrange equations. Similarly, some form of 
canonical Hamilton equations is obtained when the balance J7J is unified with 
the inverse Legendre transformation, i.e., 

f J r i i _ 

— = -g ij Pj , Of j = J[ip] 3k Y, k m g m \ 
where, obviously, J[<p] lk J[<p]kj — 5 l j (do not confuse J[ip] with g-shift of indices 

of JM). ' 

Obviously, the general balance form may be used for dissipative non-Largan- 
gian models. Simply the covariant force Qi and the generalized internal force 
Q l j (affine moment of forces, hyperforce) must involve appropriately defined 
dissipative forces (in the case of affinely-constrained continuum one can also 
consider the mutual coupling of mechanical phenomena with discrctized thermal 
effects). 

As shown in the mentioned papers, the above equations of motion may 
be formulated in various equivalent forms adapted to the kind of considered 
problems. For example, instead of the canonical (Hamilton) form, one can 
write them down in purely kinematical velocity-based terms, i.e., 

where contravariant forces F l and hyperforce N tJ (affine dynamical moment) 
may depend on all possible arguments, i.e., t, x l , tp l A, dx l /dt, dip l A/dt. Ob- 
viously, for potential models they depend only on generalized coordinates and 
possibly on the time variable t itself, and then 

/•' - ir'Qj = ■r ! '^- = <r = -^A-^-gV. (10) 

Remark: In spite of the tradition based on Riemannian geometry and relativity 
theory we shall refrain from the graphical identification of symbols F l , A™ 
respectively with Q l , Q lJ . In our treatment this would be just confusing, because 
we shall use various prescriptions for shifting the tensorial indices, i.e., various 
isomorphisms between contravariant and covariant objects. 
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As mentioned, the above equations of the motion © may be derived di- 
rectly in Newtonian terms, basing merely on the d'Alembert principle and its 
underlying spatial metric g in M. The primary quantities of this approach are 
the monopole and dipole moments of the distributions of linear kinematical mo- 
mentum and forces within the body. These quantities, just as all high-order 
multipoles may be defined for any unconstrained system of material points, 
does not matter finite, discrete, or continuous. Let $(£, a) 1 denote as previ- 
ously Cartesian coordinates of the current position of the a-th material point 
at the time instant t, x l (t) be the current position of the centre of mass, and 
!F l (t, x, a), dx/dt, (d$/dt)(t, a); a) be the density of forces per unit mass. 
As mentioned, the affine constraints are not yet assumed. 

The monopoles are simply the total quantities: the total kinematical mo- 
mentum k l (do not confuse it at this stage with the canonical one pt) and the 
total force F l affecting the centre of mass motion, i.e., 

f d& f ■ ( dx <9$ \ 

k l = J —(t,a)dfi(a), F l = J F U x, $(t, a), — , — (t, a); a) dn{a). 

The dipole moments with respect to the centre of mass current position are 
referred to as kinematical affine spin (do not confuse it at this stage with 
the canonical one £* j ) and the affine moment of forces (not to be confused 
with its potential version Q % j). They are given respectively by the following 
expressions: 

f . . I i)<\' • d V s \ 

KV = J (&(t, a) - x*) {^(t, a) - — J d/i(a), 

N** = J (^(i.a)-^)^' (t,xMt,a),^,^(t,a);a^dfx(a). 

The dipoles may be also referred to some space- fixed centre O € M, e.g., 
the origin of Cartesian coordinates in M. The difference is that "the lever arm" 
($ l - x l ) is then replaced by <P itself, and its velocity (d^/dt - dx z /dt) by 
(d&/dt)(t, a). The resulting dipoles will be denoted respectively by K(0) lj , 
N(0) l: > . For the sake of uniformity, it may be also convenient to denote the 
previous dipoles by i4T(cm) 1 -' and JV(cm)'j instead of Jf y , A™. We shall also use 
affine moments of the centre of mass characteristics with respect to the origin O. 
Thus, the translational (orbital) affine moment of kinematical linear momentum 
and translational affine moment of forces are as follows: 

dri 

K t JO) v = x l y = mx l , NtJOY 3 = x'FK 

dt 

The doubled skew-symmetric parts of the above quantities, i.e., 

S ij A"- A-. Ltr(o) ij = K tI (oy j -K tI (oy\ j{oy j =K(o) ij -K{oy\ 

tfij = N ij -N ji , Aftr(G) ij = Nt r {Oy j -Nt r {0) ji , N(0) l] = N{Oy j -N(0) ji , 



16 



represent the kinematical angular momentum and the moment of forces (torque) . 
They also occur in three versions concerning, respectively, the internal motion 
(thus S is spin), motion of the centre of mass with respect to O, and the total 
motion with respect to O. 

If now we assume that the motion is affine, then the above expressions 
simplify to 

= ^ A ^-J AB , K{Or = K tr {Or + K» =mx i ^- + ^ A ^-J AB , 
at at at 

N{0) i3 = N tl {Oy j + N ij = x i F j + N lj . 

Obviously, F and N become now functions of t, x l , dx % /dt, f z A, dip 1 A/dt. Let 
us stress that, just as it was the case with the kinetic energy, the above additive 
splitting into translational and internal parts is based on the assumption that 
the current centre of mass has permanently Lagrangian coordinates a K = 0. 
This is consistent because barycenters are invariants of affine transformations. 

By summation of elementary time rates of work over the body constituents, 
one can show that in the affine motion the total rate is given by 

fj T i 

V = g ij —Fi+g ij n i k N k i. 

Let us remind however that, besides of active generalized forces F, N controlling 
affine modes of motion, there are also hidden structural forces keeping affine 
constraints, i.e., reactions. Their density J-r does not vanish, however, their 
monopole and dipole moments Fr, Nr do because, according to the d'Alembert 
principle, the reaction time rate of work vanishes for any constraints-compatible 
virtual velocities, i.e., for any possible dx 1 /dt, Q: 

V R = 9ij-^F R i + g l3 W k N R k i = 0. 

Therefore, the effective reaction-free equations of motion are obtained from 
the primary non-constrained system by calculating the monopole and dipole 
moments. 

The above derivation is quite general and valid for all kinds of forces, includ- 
ing non-potential and dissipative ones. It relies only on the metric structure g 
in M and on the d'Alembert principle. Obviously, if equations of motion follow 
from the Lagrangian L = T — V(x, e), T given by (TTJ), then the above analysis 
implies equations (flT)]) . 

Similarly, one can easily show that 

A'-' >:',„//'"'• /■•' //''/'.• (ii) 

But these relationships become false when Lagrangian depends on velocities not 
only through the kinetic energy T but also through some generalized potential 
V, e.g., when magnetic or gyroscopic external forces are present. This is one of 
reasons we avoid denoting K 1 ^ by S y or K l j by S l j . 
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Kincmatical quantities k l , K % i are intuitive because of their direct opera- 
tional interpretation in terms of positions and velocities. At the same time they 
are lowest-order multipoles (monopoles and dipoles) of the distribution of kinc- 
matical linear momentum within the body, and it is difficult to over-estimate 
the application of multipoles and all moment quantities in practical problems of 
mechanics and field theory (cf, e.g., all Galerkin-type procedures [57]). On the 
other hand, their canonical counterparts pi, j have a very deep geometrical in- 
terpretation as Hamiltonian generators of fundamental transformation groups. 
Because of this, they are very often important constants of motions. In me- 
chanics of affinely-rigid body, equations of motion are equivalent to the balance 
laws for Pi,T, l j or, in a sense equivalently, to the ones for k l ,K^ , because La- 
grangians of non-dissipativc models, or at least Lagrangians of non-dissipative 
background dynamics, establish some link between these concepts. Similarly, 
in rigid-body mechanics equations of motion are equivalent to the balance for 
Pi , (E*,- - g ik g 3m Z m k ) or for k\ S l K 

Equations of motion ^ may be written in several mutually equivalent bal- 
ance forms. Let us quote some of them based on kinematical quantities like 
k % ,K % 3 or their co-moving representation k ,K AB , where, obviously, k l = 
ip l A_k A , K l i — (p 1 Atpi bK AB ■ The co-moving components F A ,N AB of gener- 
alized forces are given by analogous expressions, thus, F l — ip % A F A , A™ = 
^ A ^ B N AB . 

The dynamical balance expressed in terms of kinematical (non-canonical) 
quantities in spatial (Eulerian) representation reads: 

**L = F i = WaW^jab , Nij (12) 

dt dt dt dt 

For non-dissipative potential systems with Lagrangians L = T — V(x, if), it re- 
duces to ([7]), because then (TT0|) holds. Let us observe that even in the interaction- 
free case, when N = 0, the balance for K is not a conservation law due to the 
first non-dynamical term on its right-hand side. One can write that 

dt dgij 

On the Hamiltonian level, this means that the non-conservation of K even in 
geodetic motion is due to the fact that the kinetic energy depends explicitly on 
the spatial metric tensor. Affine symmetry of degrees of freedom is broken and 
reduced to the Euclidean one. 

The system (|12[) may be written in the following form: 

d ^ = F\ d J^ =m d 4 d 4 + d 4^ d -4^J AB + N { Or, (13) 
dt dt dt dt dt dt 

as a balance for the kinematical linear momentum and the total affine momen- 
tum with respect to some space- fixed origin O € M. 

If the body is rigid in the usual metrical sense, i.e., all distances between its 
constituents are constant, then the d'Alembert principle implies that the second 
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subsystems in (IT^I) . (fT5)l are to be replaced by their skew-symmetric parts, thus, 

To these equations the rigidity condition, i.e., t\ab — Qij^p 1 A 1 ^ ' b, may be auto- 
matically substituted without paying any attention to reaction forces responsible 
for the metrical rigidity. 

The above balance laws for kinematical angular momenta become conserva- 
tion laws in the interaction-free case, and even under weaker, realistic conditions 
that N or N(0) is symmetric. This is very natural on the Hamiltonian level, 
because S, J{0) are then directly related to their canonical counterparts. The 
latter are Hamiltonian generators of the isometry group of (M,g), thus, accord- 
ing to the Noether theorem, they are constants of motion in geodetic problems, 
because gyroscopic kinetic energy is isometry-invariant. 

Another very convenient balance form of equations of motion is obtained 
when generalized velocities dip 1 A/dt are expressed through the non-holonomic 
quantities then 

rlh i dK^ 

-t-vjf + w. (14) 

Expressing our balance in co-moving (material) terms we obtain 

^ = -k B J BC K CA + F A , ™£L = -K AC J CD K DB + N AB , (15) 
at dt 

or, using non-holonomic velocities, 

= — mCl A gV B + F A , jAC^Cl = _( 1 B d ^D c jCA + ^AB 

dt dt 

It is a nice feature of the co-moving representation that all non-dynamical terms 
are built only of expressions k A ,K AB or v A ,Cl A B without any direct using of 
mixed quantities like ip l A , dip 1 aI dt. The second (internal) subsystems are exactly 
afhne counterparts of gyroscopic Euler equations and exactly reduce to them 
when the rigid-body constraints arc imposed. The relationship between two 
co- moving forms is based on the equation K AB — Cl B cJ CA following directly 
from the definition of K. There is some relationship between this formula and 
Legendre transformation for Lagrangians L = T — V(x,ip). Namely, one can 
show that the internal part of |[8} may be equivalently written in the following 
form: 

t A B = k AC G CB = G BC n c D J DA , (16) 

where G £ U* ®U* denotes as previously the Green deformation tensor. There- 
fore, the canonical affinc spin is obtained from the kinematical one by the G- 
lowering of the second index. As we saw, there was a similar formula (jlip in 
the spatial representation, i.e., 

S 4 , = K lm g mj . (17) 
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It is seen that there is an easy possibility of confusion. Namely a superficial 
analogy with the last formula might suggest us to use the ^-shifting of indices 
for establishing the Legendre link between K and E. However, for any rea- 
sonable Lagrangian Y> A b ^ K ag t}cb except the special case of metrically-rigid 
motion. This is an additional reason for avoiding ambiguous symbols like T, AB 
or K b- More generally, if some tensor objects in V are related to each other 
by the g-shifting of indices, then the corresponding co-moving objects in U are 
interrelated by the G-shifting. And conversely, if two tensors in U are interre- 
lated by the 77-shift of indices, then their spatial counterparts in V are obtained 
from each other by the C-shifting, where C G V* ® V* is the Cauchy defor- 
mation tensor. The contravariant inverses of G and C are carefully denoted 
by G G U <g> U and C G V <g> V, i.e., G AC G CB = S A B , & k C k] = 5 l r The 
notation G AB , C 4 - 7 would be misleading because of the possible confusion with 
the objects r] AC r] BD GcD, g tk g^ m Ck m obtained from G, C by the usual 77- or 
^-metrical operations on indices. 

For Lagrangian systems L = T — V(x, ip) generalized forces Q, N are in- 
terrelated by (jTUJ) , thus, just as in (fT7|) . Q % j = N lm g m j. But in the co- moving 
representation, in analogy to (fl~6|) . we have Q a b = N ac Gcb, Q A b ^ N Ac ricB- 
This has to do with different (/^-transformation properties of Q, N, i.e., Q l j = 
V 1 aQ a b(v~ V ) B ^ ^ = V 1 aV 3 bN ab - similarly for E 4 3 , K*', £ a b, K AB . 

As seen from equations (|14|) . (fT5|) even in the interaction- free case neither 
K l i nor K AB are constants of motion. The same concerns their canonical coun- 
terparts E 4 j, t A B - The reason is that the kinetic energy is not invariant under 
spatial and material affine transformations (except translations, of course). At 
the same time, purely geodetic Hamiltonian models with L = T are physically 
useless because, except of rest-states, all their trajectories (straight lines in Mx 
LI([/, V)) escape to infinity. In particular, the body may expand to infinity and 
contract in finite time to a point. The metric on Q = Mx Ll(U,V) underly- 
ing the kinetic energy |T]) is unable to encode realistic interactions and predict 
elastic vibrations in purely geodetic terms. 

4 Dynamical affine invar iance 

Basing on the motivation presented in previous sections, we shall now consider 
some models which are ruled by affine groups not only on the kinematical but 
also on the dynamical level. In particular, we shall discuss left- and right- 
invariant Riemann metrics on linear and affine groups or rather, more precisely, 
on their free-action homogeneous spaces. We concentrate on geodetic models, 
when there is no potential term and the structure of interactions is encoded in 
an appropriately chosen metric tensor on the configuration space. 

Let us begin with the internal sector, when translational degrees of freedom 
are frozen and the configuration space reduces to Qint = Mx L1(U,V), or 
equivalently to F(V) (when for simplicity we put U = R"). According to the 
transformation rules for f2, f2, the most general metric tensor on Q m t invariant 
under the action of GL(V) through ([4J is that underlying the kinetic energy 
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form given by 

Tint = A D C& A B& C D , (18) 

where coefficients C are constant and symmetric in bi- indices ( B a), { D c)- This 
quadratic form is also assumed to be non-degenerate, although not necessarily 
positively definite. As £1 is a non-holonomic velocity, i.e., it is not a time deriva- 
tive of any system of generalized coordinates, the underlying metric on Qmt is 
curved. 

Quite similarly, the most general kinetic energy invariant under material 
affinc transformations ([5]) has the following form: 

Tint = tgVhrtj&u (19) 

where 1Z is also constant and symmetric in bi-indices ( 3 \), ( l k). The underlying 
metric tensor on Qj nt is also curved, i.e., essentially Riemannian. 

In general, (|18[) is not right, i.e., materially, invariant under GL(t/) acting 
through ([5]), and (fTO]) is not invariant under GL(V) acting through (g]) in the 
physical space, i.e., on the left. The exceptional situation of simultaneous spatial 
and material invariance leads us to 

A B A - ~ B - 

Tint = -Z O^jQ^i + — f2\fF j = —0, lVL k + L j 

where A, B are some constants. Using invariant terms, we can say that such 
Ti n t is a linear combination of two basic second-order Casimir invariants, i.e., 

T int = ^Tr(tt 2 ) + |(Trtt) 2 = ^Tr(£7 2 ) + |(Trf)) 2 . (20) 

Such Ti n t is never positively-definite. The reason is that the maximal semisimple 
subgroups SL(F), SL(U) (their determinants equal to unity) are non-compact, 
thus, the quadratic form Tr(f2 2 ) = Tr(f2 2 ) has the hyperbolic signature (n(n + 
l)/2 +, n(n— 1)/2 — ), where the positive contribution corresponds to the "non- 
compact" dimensions and the negative one to the "compact" dimensions in 
GL(V), GL(U). 

By the way, the above quadratic forms reduce to the Killing form (Killing 
scalar products) on L(V), L(U) [23 123 HZ| when A = 2n, B = -2. As L(V), 
L({7) are non-semisimple, in this special unhappy case the scalar product (ki- 
netic energy) is degenerate, thus, non-applicable in usual mechanical problems. 
The singularity consists of dilatational Lie algebras RIdy , Rldj/. More generally, 
the same holds when A — —Bn. Paradoxically enough, non-degenerate forms 
(|20|) (A 7^ —Bn) may be mechanically useful in spite of their non-definiteness. 

The usual d'Alembert model ([I]) invariant under additive translations (J3j> is 
the special case of general models of the following form: 

T ^~2 A 1 j ^r^T' (21) 
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where A is constant and symmetric in bi- indices ( K i), ( L j)- The peculiarity of 
(fTJ) within this class is that A factorizes, i.e., A K i L j = gijJ KL , and is invariant 
under the left action of SO(V, g) and the right action of SO(U, J) (in particular, 
SO(U, 77), when the inertia is isotropic, i.e., J = /177). It is clear that the ^4-based 
models of T; nt are never invariant under GL(V r ), GL(C/). The underlying metric 
on LI(U, V) is flat. 

Let us now consider the translational sector of motion. The only model of 
translational kinetic energy invariant under GAf(M) (afhne group of M) is as 
follows: 

„ TO „ dx % dx° m a R / \ 

T tr =— Ca— — = —Vabvv- (22) 

2 3 dt dt 2 ' y ' 

It looks like the usual kinetic energy, however, there is a very essential difference. 

Namely, in the above expression the velocity vector is not squared with the help 

of the constant and absolutely fixed spacial metric g £ V* <£> V*. Instead of it, 

the Cauchy deformation tensor C is used as an instantaneous metric tensor of 

M. Being a function of the internal configuration <p € LI([/, V), it depends on 

time. It is so as if the instantaneous internal configuration created a dynamical 

metric in an essentially amorphous affine space M. In this sense the model is 

an oversimplified toy simulation of general relativity. At the same time it is 

clear that T tr is not invariant under Gh{U) because the fixed material metric 77 

restricts the symmetry to 0(U,r)) (77-rotations of U). 

If we wish the translational kinetic energy to be GL([/)-invariant, then the 

only reasonable model is the usual one, based on the fixed M-metric g, i.e., 

to dx l dx J to ^ A ^ B 

T * = j 9ij ^r^r = j GABV v ■ (23) 

It is impossible to construct any model of T tr and T = T tt + T int , which would 
be purely affine both in M and N; in one of these spaces some metric structure 
must be assumed. Therefore, although T- lnt alone may be affine simultaneously 
in M and N, there are no reasons to stick to such models, the more so they are 
never positively-definite. These problems have to do with the non-existence of 
a doubly (left- and right-) invariant pseudo-Riemannian structure on the affine 
group GAf(n, E) ~ GL(n, M) x s M. n . Any doubly-invariant twice covariant tensor 
field on this group is degenerate. Therefore, it is reasonable to concentrate 
on kinetic energies which are affine in M and 77-metrical in N or, conversely, 
5-metrical in M and affine in N. The corresponding geodetic models in Q 
have a maximal possible symmetry, being at the same time really true geodetic 
problems (non-singular underlying metric). Such models are special cases of 
(TT5I) , p^|) . (|2"2"]) . (|2"5|) . thus, we start with some statements concerning the general 
case. 

For the model (p~8|) . (|22|) affine in space and matrical in the material, Legen- 
dre transformation has the form: 

t A B =£ A b C d£i D c, Pa — ""Iab^i (24) 
where, obviously, the second equation may be rewritten as 

m( ', ,! '. (25) 
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The corresponding geodetic Hamiltonian is given by T = 7^ r + 7j nt , where 

T tr = ^ v AB p A p B = ^CVpiPj, r tat = ^s^E^E^c, 
2m 2m 2 

and the symmetric bimatrix C is reciprocal to £, i.e., 

rA K r L C _ sA rC 
i- B LJ~ K D — DO B- 

For the model IT51) . (f2"3")l . metrical in space and affine in the material, Leg- 
endre transformation may be represented as follows: 

E^ftVifi'*. Pi = m gij v j , (26) 

where, in analogy to the £-case, the second subsystem may be rewritten as 

Pa = 7nG AB v B . 

Inverting the Legendre transformation, we obtain geodetic Hamiltonian T = 
%,t + Tint, where, dually to the £-formulas, 

Ttr = = ^ ABpApB ' %nt = l&Si&i&k. 

Obviously again 1Z denotes the inverse bimatrix of K, i.e., K a h k {R} k h = 8 a l 8\. 

Let us admit non-geodetic models of the form L = T — V(x,ip), H = 
T+V(x 1 ip), where V(x, if) is a usual potential energy depending only on the in- 
dicated configuration variables. Then the balance equations for £-models (affine 
in space, metrical in the material) read: 

^ = Q„ ^-=--C ik p kPj +Q i j , (27) 
at at m 

with the same meaning of generalized forces as previously, i.e., 

dV dV 

When taken together with the Legendre transformation, the above balance laws 
are equivalent to the Hamilton canonical equations. They may be also gener- 
alized so as to include some non-Hamiltonian, e.g., dissipative terms on their 
right-hand sides. It is seen that in the purely geodetic case (when Qi = 0, 
Q l j = 0) the canonical linear momentum is conserved, but the affine spin is 
not so due to the first term on the right-hand side of the balance law for E. 
The reason is that S J j generate linear transformations of internal degrees of 
freedom; these transformations do not affect translational variables. Therefore, 
the full affine symmetry is broken, and E is not a constant of motion. But 
one can reformulate the balance laws (|2"T|) , just as in the d'Alembert model, by 
introducing the total canonical affine momentum 1(0), related to some fixed 
origin O £ M, 

nor, : A(oy 3 + e^ = x*pj + sy 
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The first term refers to the translational motion, the second one to the relative 
(internal) motion. Something similar may be done for generalized forces, 

Q„,(OY, : 0tr(O)*i + Q l j = ■>■'(■>, + Q i j- 

Then the system of balance equations l|27p may be written in the following 
equivalent form: 

It is seen that in the geodetic case one obtains conservation laws for pi, l(0) l j, 
i.e., for the system of generators of the spatial afhne group GAf(M), just as 
expected. 

Let us observe some funny feature of our geodetic equations. Namely, the 
canonical linear momentum is a constant of motion, but the translational ve- 
locity is not. This is because the Legendre transformation (1251) implies that the 
translational motion is influenced by internal phenomena. Except some spe- 
cial solutions even the direction of translational velocity is non-constant and 
depends on what happens with " internal" degrees of freedom. This is a kind of 
"drunk missile" effect. Something similar occurs in the dynamics of defects in 
solids 28J. It is also non-excluded that the non-conservation of velocity might 
be an over-simplified model of certain specially-relativistic phenomena (internal 
motion results in changes of internal energy, and therefore, in the rest mass 
pulsations; but the latter ones influence the effective inertia, and therefore, the 
translational motion). One can show that the time-rate of translational velocity 
may be expressed as follows: 

m ^t~ = -c aj ^v b + F a = mv b (Q\ + c ad n m d c mb ) + F a , 

where the contravariant force F is given by the following expression: F a = 
C ab Qb (because of the formerly mentioned reasons, we avoid denoting F a as 
Q a ). It is explicitly seen that v is variable even in the purely geodetic motion. 
The £-based geodetics in Mx LI(/7, V) do not project onto straight lines in M. 

Let us now consider the 7?.-based models (fT9]l. ([23]). metrical in space and 
affine in the material. As mentioned, they are somehow related to the Arnold- 
Ebin-Marsden-Binz approach to the dynamics of ideal fluids [H HI 02 [13l HH [331 
[34] . Our Poisson brackets imply that the balance form of equations of motion 
may be expressed as follows: 



dp a 

again with 



Q a = -g, Q\ = -^t/b = (v-YiQV* 

in the potential case. In geodetic models p a , Ti A b are conserved quantities as 
explicitly seen from the balance equations and expected on the basis of invari- 
ance properties. Indeed, the 7?.-model of T is invariant under the Abelian group 
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of translations in M. Therefore, p a are constants of motion as Hamiltonian 
generators of this group. Similarly, as seen from our Poisson brackets, S A s 
Poisson-commute with p a and S 1 ^, therefore, with the total geodetic Hamilto- 
nian T. This is due to the explicitly obvious invariance of the 7?.-based T under 
the group of material linear transformations GL([/). Surprisingly enough, the 
co-moving components of linear momentum, pa = Pi<p l a, are not constants of 
motion. But this is also clear because the material space N has in our model a 
distinguished point, i.e., the Lagrangian position of the centre of mass. Because 
of this, translations in N fail to be symmetries and their Hamiltonian generators 
Pa are non-conserved. According to the structure of Legendre transformation, 
translational velocity v a — dx a /dt — g ab pb is also a constant of motion, just as 
p itself. Therefore, in geodetic 7£-models, the geodesic curves in Mx Ll(U, V) 
project to M onto straight lines swept with constant velocities (uniform mo- 
tions). This means that there is no "drunk missile effect" and contravariant 
representation of the translational balance takes on the usual form: 

As previously, the balance laws ([29]) become a closed system of equations 
of motion when considered jointly with the Legendre transformation (|26p . Let 
us observe that the structure of (|29|) is in a sense less " aesthetical" than that 
of (|27|> because it is more non-homogeneous. The point is that in (|27p both 
subsystems are written in terms of spatial objects, whereas in (f2"9"| one uses 
the mixed representation: spatial for the translational motion and material 
for the internal one. Obviously, ([29]) may be done symmetric, dual to (|27|) . by 
substituting^ = piip 1 a- But this immediately makes the translational equation 
more complicated. 

There is also another problem. The simplicity of our balance laws (conserva- 
tion laws in the geodetic case) is rather illusory. The point is that, as mentioned 
above, the total system of equations of motion consists of the balance laws and 
Legendre transformation. The balance (f27l) . (|28|) looks simple in Euler (spa- 
tial) representation, but the corresponding Legendre transformation is simple 
in Lagrangian (material) representation (|24p . And quite conversely, the internal 
part of (|29p is simple in the co-moving terms, but its Legendre transformation 
is simple when expressed in the spatial (Eulerian) form (|26|) . 

One can easily show that the internal parts of Legendre transformations 
(f24|) . (j26|) may be respectively expressed as follows: 

= Cj k in l k , t A B = n A B c D h D c, (30) 

where 

n A B C D = (<p- 1 ) A i<P i B(<p- 1 ) C k<(tDK i j k l. 

Obviously, this form is rather complicated because the coefficients at and Cl 
are non-constant; they depend on the internal configuration ip. Simplicity of 
the balance laws is incompatible with simplicity of Legendre transformations. 



25 



As mentioned, when translational degrees of freedom are taken into account, 
there are no sensible models which would be affine simultaneously in space and 
in the material. The highest symmetry of mathematical interest and at the 
same time physically reasonable is that affine in space and rotational in the 
material, and conversely, Euclidean in space and (centro-)affine in the body. 
The latter model is an over-simplified discretization of dynamical systems on 
diffeomorphisms group as used in hydrodynamics and elasticity. 

In materially isotropic C- models the quantity C a b C d is a linear combi- 
nation of tensors rj Ac rjBD, 5 A d$ C b, S A b5 c d- Similarly, in spatially isotropic 
7?.-models the tensor TV k i is a linear combination of terms g zk gji, S l iS k j, S l jS k i. 
Therefore, (TT5)) . (fTi?|) take on, respectively, the following forms: 

^nt = ^mLn K M Ci L NV MN + ^n K L n L K + ^n K K n L Ll (31) 



T lnt = V I t}<->',<>', • f <->'.!!',• (32) 



where the constants /, A, B are generalized internal inertia scalars. It is clear 
that if I — 0, then these expressions become identical. The /-terms break 
the centro-affine symmetry in U and V, and restrict it to the metrical one, 
respectively, in the sense of metric tensors rj or g. The first term in (|31[) . just 
as (|22[) , may be expressed in terms of the Cauchy deformation tensor, i.e., 

2 3 dt dt ' v ' 

Let us observe that the isotropic inertial tensor lr\ AB in ([33]) might be re- 
placed by the general one, 



1 <kp l A dip> B jAB 1 Aif AL jAB ( n A , 

2° ij dt dt J -2 r]KLQAQBJ ■ (34) 



This expression is structurally similar to the d'Alembert formula (fT]). The differ- 
ence is that the fixed metric g is replaced by the (^-dependent Cauchy tensor C. 
There is not only formal similarity but also some asymptotic correspondence be- 
tween (fMf and |T]). Obviously, for the general J, (|3~i|) is not metrically isotropic 
and its internal symmetry is reduced to 0(U, rj) PI 0(U, J). The /-terms in (pH]) , 
(|3"2"|) are positively definite if / > 0. Moreover, the total expressions (|3"Tj) . (|32p 
are positively definite for some open range of (/, A, B) el 3 . Roughly speaking, 
the absolute values of A, B must be "sufficiently small" in comparison with /. 

The internal part of Legendre transformation (f2"l)l for £-models becomes 
now (i.e., for (f3"T| ) 

t K L = It 1 km 11ln VL N m + AVL K L + B5 K L Cl M M- (35) 

This may be alternatively written as follows: S l j = IC lb Cj a fl a b + ASl % j + 
B5 l j£l m m , it is the very special case of ([50]). 
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The inverse Legendre transformation has the same structure, i.e., 

it L = —T) T]ln2-> M + — ±> L + —0 M, (ob) 

I A B 

where I = (I 2 - A 2 ) /I, A= (A 2 - I 2 ) /A, B = - (I + A) (I + A + nB) /B. 
When written in Eulerian (spatial) terms, this formula becomes 

Si — — O Ljj a Zj {, -| — Zj j -\ —0 Zj m , 

/ A B 

with the same as previously meaning of modified inertial coefficients /, A, B. 

Similarly, for 7?.-models based on (j3"2")l the internal sector of Legendre trans- 
formation has the following form 

= Ig im g jn n n m + Arfj + BS^Q" 1 ^ (37) 

which is inverted as 

Q l j = —g lm gj ri Y, n m + -—E'j + —5 l jTj rn m . (38) 
/ A B 

The co-moving representation of these formulas is given by the following expres- 
sions: 

t K L = IG KM G LN Cl N M + ACl K L + BS K L Cl M M , 

" L — ~<J (JLN^ M + — ii L + —0 M- 
J AS 



The general balance laws (f27|). ([29)1 considered jointly with these Legendre trans- 
formations (including the obvious translational sector) provide the complete 
system of equations of motion (naturally, the definitions of si, fi are to be sub- 
stituted). These equations, due to the very special structure of (f3Tj) . ([32]) are 
relatively readible and effective. At the same time one can show that for incom- 
pressible bodies even in the purely geodetic case (Qi = 0, Q l j = 0) there exists 
an open set of solutions which are bounded in the internal configuration space 
U(U,V), so the elastic vibrations may be encoded in the very kinetic energy 
(Riemann structure) without the explicite use of forces. 

After substituting the above inverse of Legendre transformations to kinetic 
energies (|3"T|) . (f3"2"j) . we obtain the following geodetic Hamiltonians of internal 
motion: 

Tint = -Lvkl£ K m£ L nV MN + ^t K L t L K + ^Lt K K t L L , (39) 
21 2A 2B 

7I n t = iflifcE^EV + -^f&j&i + -i^EV (40) 
21 2A 2B 

In certain problems it may be convenient to write down the first formula in spa- 
tial terms; similarly, the second one may be expressed with the use of co-moving 
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representation. Therefore, we obtain, respectively, the following expressions: 



int — — C-feiZj m L „G H -2j ;2j fc H — 2j fclj z, 

2/ 2A 25 

Ant — — '■JKL^ N 1 ^ H AT H ~ ^ L- 

21 2A 2B 

The corresponding velocity-based formulas (f3Tj) , |32|) for kinetic energy may be 
written in an analogous way. Simply 1/1, 1/A, 1/B in the last expressions are 
to be replaced by /, A, B, and simultaneously one must substitute f2 fe z, Cl K l 
instead of Y. k h t K L . 

For certain purposes it is convenient to rewrite geodetic Hamiltonians (|39|) . 
(|40p in an alternative form: 

%* = ^ Tr ^ 2 ) + ^( Trf: ) 2 -^ Tr ( y2 )' W 

%nt = ^Tr(S 2 ) + i-(TrI]) 2 -^Tr(^), (42) 
la Zp 4/i 

where a := I + A, := -(I + A) (I + A + nB)/B, fx := (I 2 - A 2 )/ 1, and V, 
S denote, respectively, the vorticity and spin given by ©■ It is clear that the 
only distinction between expressions (|4ip and (|4"2")) is that in their third terms, 
thus, we can rewrite them concisely like 

%* = ic(2) + ^C(ir-^Tr(V% (43) 

T * - h c ^ + h c{1? ~h^ (44) 

where C(k) denotes the fc-th order Casimir expression built of generators, i.e., 
C(Jfc) := Tr(S fe ) = Tr(£ fc ). 

By analogy with the physical 3-dimensional case the quantity — Tr(S* 2 ) may 
be interpreted as the doubled squared norm of the internal angular momentum, 
-Tr(S* 2 ) = 2\\S\\ 2 . Similarly, -Tr{V 2 ) = 2\\V\\ 2 . Therefore, the formulas (jig]). 
(|44[) may be respectively written in the following intuitive and suggestive way: 

%* = ic(2 ) + ^C(lf + l-\\Vr, (45) 

%* = ^C(2) + ^C(1)' + £||S||». (46) 

Obviously, for the model (|20[) invariant under the left and right action of linear 
groups we have 

T ~ = ^ + m^Am C{1)K (47) 

When performing computations, it is convenient to use orthogonal coordi- 
nates, rjKL = $kl, 9ij — Sij, and rewrite some of the above formulas in terms 
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of purely matrix operations. Thus, (pTLj) . ([52"]) become, respectively, 

T int = ^Tr(M) + |Tr(r! 2 ) + |(Tr^) 2 , (48) 

T in t = ^Tr(^) + |Tr(ft 2 ) + |(Trft) 2 . (49) 

Obviously, the second and third terms in these formulas are pairwise equal 
because Tr(ft 2 ) = Tr(f2 2 ) and Tr(Q) = Tr(fi). 

Legendre transformations (f35|) . (|37| and their inverses (|36)) , p8|) are as fol- 
lows: 

S = IfF + Af2 + B(TrQ)I„, £ = 7fi T + AO + B(Trfi)I n , 

= is T + — £ + i(Tr£)/ n , n = is T + is + i(Tr£)I n , 
IAS 7 A B 

where I n denotes the n-th order identity matrix. 

Similarly, for the kinetic Hamiltonians (|39|) . ([40]) we have, respectively, 

T hlt = ^Tr(£ T £) + ^Tr(£ 2 ) + ^(Tr£) 2 , (50) 
21 2A 2B 

T int = ^Tr(E T E) + ^Tr(E 2 ) + ^(TrI]) 2 . (51) 
21 2A 2B 

This matrix representation is very lucid and useful in calculations. Never- 
theless, in comparison with the systematic tensor language, it may be risky and 
misleading because it obscures the geometric meaning of symbols and concepts. 
And this is worse than the lack of aesthetics; when no care is taken, this may 
lead simply to logical and numerical mistakes. 

We finish this section with some geometric remarks. 

Kinetic energy T of a non-relativistic mechanical system is equivalent to some 
Riemannian structure Y on its configuration space Q. In terms of generalized 
coordinates and velocities, we have that 

T = ¥ ap{q) ^^' r = r a p(q)dq a ®d q e. 

Although usually somehow related to the metric tensor g of the physical space 
M, Y need not be directly interpretable in terms of geometrical distances in M. 
As a rule, it depends not only on g but also on certain parameters characterizing 
inertial properties of the system, i.e., masses, inertial moments, etc. 

It is instructive to describe explicitly in a bit more geometric form the metric 
tensors Y on Q — M x Ll(U, V) underlying the kinetic energies defined above. 
For this purpose it is convenient to introduce auxiliary geometric objects. 

Let Ea, &% denote the basic vectors in U, V underlying our Lagrange and Eu- 
ler coordinates a K , xK The corresponding dual basic covectors in U* , V* will be 
denoted as usual by E A , e\ Generalized coordinates in Q will be, as previously, 
denoted by x l , a, and no misunderstandings just simplifications follow from 
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using the same symbol x l for coordinates in M and for their pull-backs to Q. 
Now, we introduce two families of Pfaff forms (differential one-forms) on Q, i.e., 
uj A b '■= ((f~ 1 ) A idip l B, u l j :— d(p t A{ ( p~ 1 ) A j- These basic systems depend on the 
choice of bases E, e, but this is, so to speak, a covariant dependence. In other 
words the L([/)- and L(V)-valued one-forms Cj := lj a bEa®E b , cj :— o/j-et®^ 
are base-independent. 

In addition, we shall need the following two families of Pfaff forms: 6 — 
{ip~ 1 ) A i dx l , 9 l — ip l a8 A = dx l . Just as previously, they give rise to the objective 
base-independent U- and V- valued differential one- forms: 9 — 9 a Ea, 9 = 9 l e,i. 

The base-independent objects u),6 and u>,9 could be defined without any 
use of bases, however, the above definitions are technically simplest. The above 
objects are closely related to the concept of affine velocity in co-moving (La- 
grange) and spatial (Euler) representations. Namely, for any history R 3 t h 
(x(t) , ip(t)) , the quantities £1 A b, are evaluations of u) A b> ^ l j on the tan- 
gent vectors (general velocities) given analytically by dx l /dt, d<p l A/dt. Roughly 
speaking, we could say that 

This is obviously a kind of joke, but fully justified on the basis of infinitesimal 
Leibniz notation. Similarly, the co-moving and spatial components of transla- 
tional velocity are given by evaluations of 9 A and 9 l on tangent vectors, and 
using the same trick we could say that 

v = — , V 1 = —. 
dt dt 

At any point of the configuration space the systems 9 t ,oj k i and 9 a ,cu k l 
provide two bases in the corresponding space of covariant vectors. We could 
ask for the corresponding bases of contravariant vector spaces. It is convenient 
to use the language of contemporary differential geometry, where vector fields 
X with components X % (meant in the sense of some local coordinates q 1 ) are 
identified with first-order differential operators of directional derivatives Vjc, 
i.e., X — X % d/dq % . One can easily show that the bases Hri E A b and Hk, E a i,, 
dual respectively to 9 K ,(2> a b and 9 k ,u> a b, are given by the following differential 
operators: 

TJ — ,j ^ TT — P A - ,n k d ™ _ a 9 

OX OX K Of A Oip°K 

The general C- models (TT51) , (1^1) are based on metric tensors of the following 
form: 

T = mriAB6 A ® 9 B + C b a d c^ A b ® tb°D- 
Similarly, for 7\L-models (fT§)) . we have 

L = mg i j6 i ® 9 3 + W ® u k i- 
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The corresponding contravariant metrics underlying the kinetic Hamiltonians 
are given, respectively, for £-models by 

f = - V AB H A ®H B + C B A D C E A B ® E C D 
m 



and for 7?.-models by 



1 



r = -g^Hi ® H 3 + TVh&j ® E k - 



m 



If the kinetic energy of internal motion (|2"0)) is invariant simultaneously under 
GL(V / ) and GL(£7), then the corresponding metric tensor on Ll(U, V) is given 

by 

rf nt = Auj k l ® lo L k + BCj k k <g> Cj l l = Auj k i ® J k + Bw fc fc ® w\ 
and its inverse by 

= J Ek ^ Elk -A(A^B) Ek ^ El - 

Here rP nt is a linear combination of two Casimir-like objects built of Pfaff forms 
w in a quadratic way. As already mentioned, r? nt becomes the group-theoretic 
Killing tensor when A = 2n, B = —2. This is just the pathological situation to 
be excluded because for A/B = —n the tensor T° nt is singular. 

For our models afhne in space and metrical in the body we have that 

r = m VKL e K <g> § L + i mL r, MN u K M ® u L N + r° nt . 

Similarly, for models metrical in space and affinc in the body: 

r = m giJ e l ® e j + i<j,-,</ , ® u k i + rf nt . 

The corresponding contravariant (reciprocal) metrics are given by 

f = -r] KL H K ®H L + hn KL r] MN E K M ®E L n +Le k L ®E L K + Le k k ®E l l 
mi A B 

for spatially afhne models and 

f = —g**Hi O Hj + —//,/.. '/'''/. ' , 8) £ fc z + ® + ® £^ 

m J A B 

for materially afhne models. Obviously, the last two terms in both expressions 
coincide. The kinetic (geodetic) terms of Hamiltonians for geodetic and poten- 
tial systems are based on T-tensors, namely, 

T = ~T^(q) Pfl p v , 
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where r MK r Kt/ = and p M = dL/dq^ = dT/dq^ 1 are canonical momenta 
conjugate to q^ (dual objects to generalized velocities <p). 

The above objects to, 9 and B, 7? possess natural generalizations to curved 
manifolds with affine connection. They appear there, respectively, as the con- 
nection form, canonical form, fundamental vector fields, and standard horizontal 
vector fields on FM (the principle fibre bundle of linear frames in a manifold M) 
[27] . Such a formalism is used in mechanics of infinitesimal affinely-rigid bodies, 
when affine degrees of freedom are considered as internal ones, i.e., attached to 
material points moving in manifolds with curvature and torsion [ST]. 



5 Without translational motion 

It is instructive to consider the simplified situation when the centre of mass is 
at rest and the covariant translational forces do vanish, i.e, Qi = 0. For the 
"usual" d'Alembert model such a situation, characteristic for practical elastic 
problems, was discussed detailly in our earlier papers [151 H51 H51 [521 1531 I5H 155] . 
And from the purely geometric symmetry point of view nothing particularly 
interesting happened there due to this simplification. In our model, based on 
dynamical affine symmetries, the translationless situation is an important step 
of the general analysis. 

We have mentioned that in spatially affine £-models with Qi — 0, in par- 
ticular in geodetic ones, canonical linear momentum is a constant of motion 
but translational velocity is not (except some very special solutions). This vi- 
olates our ideas about Galileian symmetry, at least in the form developed in 
the "usual" d'Alembert mechanics. Nevertheless, the concept of translationless 
motion is well-defined because in the usual potential models equations v l = 0, 
Pi = are equivalent; this is one of exceptional cases when the constancy of 
velocity does not contradict the constancy of linear momentum. One must only 
remember that the Galilei transforms (in the usual sense) of such space-resting 
solutions will not be solutions any longer. 

In £-modcls without translational motion the evolution is ruled by the second 
of the balance laws ([2"T]) with the simplified right-hand side, i.e., 

dt ^ r 

Affine invariance in M implies that in the completely geodetic case this becomes 
simply the Noether conservation law: 



i.e., affine spin in spatial representation is a constant of motion. As mentioned, 
to obtained a closed system of equations, one must consider the above balance 
(conservation) jointly with the Legendre transformation and the definition of 
affine velocity. Unfortunately, the nice form ([24]) with constant coefficients 
cannot be used because T* A b is not a constant of motion in the geodetic case. 
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And in general, (|50|) is too complicated to be effectively used. But it turns 
out that something may be done for our simplified model (|3ip . affine in M and 
^-metrical in U. 

Something similar may be said about 7?.-models, moreover, they are in some 
respects simpler. The balance equations (|29[) reduce to their internal parts, i.e., 

~dT = Q B ' 

and become conservation laws for the co-moving affine spin in the geodetic case, 
i.e., 

dt A R 

= 0. 



dt 

As previously, the simplicity is only seeming one because the laboratory compo- 
nents T, l j fail to be constants of motion and Legendre transformation expressed 
in co-moving terms (|30[) . in general, is rather complicated. Fortunately, for our 
models (|32[) , metrical in space and centro-affine in the material, also something 
may be done. 

Let us begin with the over-simplified model with I — 0, affine both in 
the spatial and material sense. It is easily seen that the general solution for 
translation-free geodetic motion is then given by the system of orbits of one- 
parameter subgroups of GL(V) or, equivalently, one-parameter subgroups of 
GL(U), i.e., 

(p{t) = exp(Et)(p = ipo exp(Et), (52) 

where ipo is an arbitrary element of Ll(U,V), E is an arbitrary element of 
L(V) = GL(V)', and E = (p' a 1 E(p is the corresponding element of L(U) — 
GL(t/)' obtained by the <^o-similarity. If we identify formally U and V with R" 
(by the particular choice of bases), then the phase portrait consists of all one- 
parameter subgroups of GL(n,R) and of all their left cosets or, equivalently, 
of all their right cosets. One must only remember that although the sets of 
left and right cosets coincide, they are parameterized in a different way by 
the corresponding generators and initial shifting elements. The reason is that 
GL(n,R) is non-Abelian and, in general, its one-parameter subgroups are not 
normal. If we write the group-theoretical version of (|52p . i.e., 



g(t) = exp(at)h — hexp(h 1 aht), 

it is seen that the coinciding left and right cosets usually refer to different 
generators a and h~ 1 ah, thus, to different subgroups. If some left and right 
cosets refer to the same subgroup, i.e., the same generator a, and have non- 
empty intersection, then, as a rule, they are different subsets, i.e., 

gi(t) — exp(at)h ^ /iexp(at) = t/2(i)- 

Only the dilatational subgroup is exceptional because, consisting of central ele- 
ments, it is a normal subgroup, and h~ 1 ah = a for any dilatation generator a; 
the previous inequality becomes equality for any h G GL(n, R). 
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Let us notice that in the solution ([52")) the pairs (fo,E and ipo , E play the 
role of differently represented initial conditions; in this sense they label the 
general solution. Thus, (po = <f(0) is an initial configuration, whereas E = 0(0), 
E = 0(0) are initial and at the same time permanently constant values of 
the laboratory and co-moving afhne velocities. Therefore, the initial values of 
generalized velocities are given by <p(0) = Eip^ = tpoE. 

It is seen that for 1 = the structure of general solution resembles that 
of the spherical rigid body. It is so for every geodetic model on a semisimple 
group or its trivial central extension if the kinetic energy is doubly (left and 
right) invariant 2J. But we should remember that even in the simple case of 
a free anisotropic rigid body situation changes drastically. Kinetic energy is 
invariant under left regular translations on SO(3,R) (identified with the config- 
uration space) but no longer under right translations. As a rule, one-parameter 
subgroups and their cosets fail to be solutions, i.e., they are not geodetics of 
left-invariant metric tensors on SO(3,R). There are some exceptions, however, 
namely stationary rotations 2 , 35 , 36 , 42J . They happen when one of main axes 
of inertia has a fixed orientation in space and the remaining two perform a uni- 
form rotation about it with a fixed angular velocity. Thus, there is a subset of 
general solution given by three one-parameter subgroups and all their left cosets 
(the non-moving axis of inertia may be arbitrarily oriented in space) . This is the 
special case of what is known as relative equilibria [TJ [351 G23 • They correspond 
to critical points of geodetic Hamiltonians restricted to co-adjoint orbits in the 
dual space of the Lie algebra SO(3,R)' ~ SO(3,R) [DOS [36]. Such particular 
solutions, although do not exhaust the phase portrait, contain an important 
information about its structure. 

Something similar happens in our afhne models when I =/= 0. The general 
solution is not any longer given by orbits of one-parameter subgroups. Neverthe- 
less, there exist geometrically interesting orbits which are particular solutions, 
i.e., generalized equilibria. 

Let us begin with geodetic £-models ([3~Tj) left-invariant under GL(V) and 
right-invariant under 0(U,rj). One can show after some calculations that there 
exist solutions of the following form: 

<p(t) = ^oexp(Ft), (53) 

where the initial configuration (po G IA(U,V) is arbitrary just as in f|52[) . But 
now F € L(J7) ~ GL({/)' is not arbitrary any longer, instead it must be 77- 
normal in the sense of commuting with its ^-transpose, i.e., F A cT] CD F e dt\eb — 
r\ AD F E dVecE c b = 0. Introducing the 77-transpose symbol, 

(F^ t ) a b :=V AC F D cVdb, (54) 

we can write the above condition in the following concise form: 

[F, F" T ] = FF nT - F r ' T F = 0. (55) 

It is obvious that for such solutions affine velocities are constant and given by 

n = (p Fip„ 1 , Ci = f. 
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The initial data ipo, F are independent of each other. The only restriction is that 
of 77-normality imposed on F alone. This holds, in particular, in two extremely 
opposite special cases when F is 77-skew-symmetric or 77-symmetric, i.e., 

pvT _ _p pvT _ p 

In the skew-symmetric case the one-parameter group generated by F consists 
of ^-rotations, i.e., 

exp(Ft) e SO(U,rj) C GL(U). 

If F is 77-symmetric, then so are transformations exp(Ft); they describe pure 
deformations in U in the sense of 77-polar decomposition. 

In calculations one identifies usually U and V with R" and their metrics 77, 
g with the Kronecker delta. Then the solutions ([53")) become all possible left 
cosets of one-parameter subgroups of GL(n, R) generated by all possible normal 
matrices F £ L(n,R), i.e., [F, F T ] — FF T — F T F — (in this formula we mean 
the usual matrix transposition) . 

Following (|52|) we can try to rewrite (|53|) in terms of the left-acting one- 
parameter subgroups. It is easy to see that 

<f(t) = <Po exp(Fi) = exp(Ft)(f , (56) 

where F = ipoF^ 1 £ L(V) = GL(V)'. 

In this representation ipg is still arbitrary but F is subject to some restrictions 
following from (|55| and depending on i^o- Namely, F is normal in the sense of 
the Cauchy deformation tensor C[<£>o] used as a kind of metric in V, i.e., 

F\C[cp } ak F l k C[<po}ij - C[cp } lk F l k C[<po}iaF a 3 - 0. 

Introducing in analogy to (f54|) the C[y?o] -transpose of F, i.e., 

(pCl^Ty. ._ c[y Q f k F l k C[y ] lh 

we can write simply 

[f,F c H r j=0. (57) 

Therefore, in the right-cosets representation the initial configuration ip Q and 
the generator F are mutually interrelated. Namely, if ipa is not subject to any 
restrictions, then F satisfies the condition ([57]) explicitly depending on ipo. And 
conversely, if F is arbitrary, then the initial conditions of ipo must be so suited to 
any particular choice of F that the commutator condition (|57j) is non-violated. 

Let us observe that in all £-models the spatial metric g does not occur in 
expressions for the kinetic energy at all. Thus, as a matter of fact it does not 
need to exist at all and the physical space M may be purely affine. Only the 
material metric rj in the body is essential for (|3~lj) . Let us notice, however, that if 
both g £ V*®V* and r\ S U*®XJ* are fixed, then some family of special solutions 
may be distinguished, for which the relationship between initial configurations 
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and infinitesimal generators is simpler and more symmetric. Namely, we can 
start from the very beginning with the representation 

<p(t) = exp(Et)(p , 

where E and ipo are respectively some elements of L(V) and Ll(U, V). It is easy 
to see that, when some metric g is fixed in V, then we have at disposal a very 
natural family of solutions assuming that ipa £ LI(C7, V) is an isometry and E 
is g- normal, 

[E, E 9T ] = EE gT - E gT E = 0, (58) 
where in a full analogy to the previous expression we use the definition 

Obviously, such solutions form a submanifold of the family (|55|) , ([57]) because 
then C[(f ] = g- 

Now let us consider geodetic TZ- models (|3^|) . which are left-invariant only 
under 0(V, g) but right-invariant under the total GL(U). Now, as expected, the 
situation will be reversed. Let us assume solutions in the right-coset form: 

ip(t) = exp(Et)(p , 

where (fo £ Ll(U,V), E £ L(T^). It is easy to show that such a curve (right 
coset) satisfies, in fact, equations of geodetic motion if E is (/-normal, just as in 
(|58[) . but now tpo may be quite arbitrary isomorphism of U onto V. And if we 
write the above curve as a left coset, i.e., 

(p(t) = ip exp(Et), E = ipo l E(p e L(U), 

then it is easy to see that, with a still arbitrary tp , E will be G[</?o] -normal in 
the sense of Green deformation tensor G[yo] G U* ® U* , i.e., 

[E, E G[ipo]T } = E E G[ipo]T - E G[vo]T E = 0, 

where the G[<£>o] -transpose is defined in a full analogy to the above r\- and g- 
transposes, 

{e gmt )Ab G[<po] AC E D c G[<p ] DB . 

Just as previously, we can distinguish an interesting submanifold of such 
solutions when some material metric tensor 77 S U* <S> U* is fixed (we know it 
does not exist in ([221) )■ They are given by curves of the following form: 

ip(t) = ip exp(Ft), 

where (po € 0(U,rj; V,g) is an isometry and F £ L(C/) is 77-normal in the sense 
of (|54p . ([55]) . For such solutions we have G[(^o] = V- Manipulating with 77 we 
introduce some kind of parametrization, ordering in our manifold of relative 
equilibria. 
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The above particular solutions are very special, nevertheless very impor- 
tant. Their position in our model is analogous to that of stationary rotations in 
rigid body mechanics. They provide a kind of skeleton for the general solution. 
Nevertheless, some, at least qualitative, rough knowledge of the phase portrait 
would be mostly welcome. The crucial question is to what extent the purely 
geodetic models may predict bounded motions. Obviously, this is impossible 
for compressible bodies, when the configuration space is identical with the total 
LI(f7, V"). To see this it is sufficient to consider the special case n = 1, when 
compressibility is the only degree of freedom of internal motion. There is only 
one affinely-invariant model of T; nt . The resulting trivial geodetic model pre- 
dicts, depending on the sign of the initial internal velocity, cither the infinite 
expansion or contraction, although in the latter case the object shrinks to a 
single point after infinite time. The only bounded (and non-stable) solution is 
the rest state. Something similar occurs in n-dimensional geodetic problems. 
Namely, degrees of freedom of the isochoric motion are orthogonal to the pure 
dilatations and completely independent of them. Some purely geometric com- 
ments are necessary here. Namely, if N and M are purely amorphous affinc 
spaces, in particular no metrics rj, g are fixed in U, V, then their volume mea- 
sures are defined only up to multiplicative constant factors. They are Lebesgue 
measures, i.e., special cases of Haar measures invariant under additive Abelian 
translations in U, V (in N, M). Let us denote some particular choices respec- 
tively by vjj, vy . Obviously, for any measurable domain Y C U and for any 
configuration tp £ LI([/, V) we have 

M<P(Y)) = AM^(F). 

The scalar multiplicator A(tp) depends on tp and on non-correlated normaliza- 
tions of vy, vjj but does not depend on Y. Obviously, for any A £ GL(J7), B £ 
GL(y) we have 

A(AipB) = (det A)A(ip) detB. 

The motion is isochoric if A is constant in the course of evolution. Obviously, 
this definition is independent of particular normalizations of vy , vjj. The mani- 
fold L1(U, V) becomes then foliated by (n 2 — l)-dimcnsional leaves consisting of 
mutually non-compressed configurations. Every such leaf establishes holonomic 
constraints, and the total foliation is what is sometimes referred to as semi- 
holonomic or quasi-holonomic constraints. If some metric tensors r\ £ U* (&U* , 
g £ V* ® V* arc fixed, then the measures vu, vy may be fixed respectively as 
v v , i/ g , and in terms of coordinates 



dv n — \/det [rjK^da 1 ■ ■ ■ da n , dv g — \J det^Jdx 1 ■ 



■dx r - 



Using Euclidean coordinates we can simply put A (ip) — det^ 1 ^-] but, obvi- 
ously, this convention fails for general coordinates. For non-Euclidean affinc 
coordinates we have 

a / \ y/det[g a ] . 
sj det \r\AB \ 
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Let us remind that the corresponding curvilinear formula reads [15) 



dv g (x{a)) = y/det[g k i(x(a))} 
du v (a) ^dct [rj A B{a)\ 



da K 



(59) 



If some volumes are fixed in U and V, e.g., due to some choices of metrics 77, 
g, then the volume extension ratio A(ip) is uniquely fixed. In certain formulas 
it may be convenient to use the additive parameter a((p) instead of the multi- 
plicative one, i.e., A((p) = exp[a(ip)]. Another convenient dilatation measures 
are those describing the linear size extension ratio, 



D(<p) = ^P) = 



exp 



a(tp) 



exp [q(ip)] 



The only possibility of stabilizing dilatations is to include some extra po- 
tential preventing the unlimited expansion to the infinite size and asymptotic 
contraction to the point-like object. There is plenty of such phenomenological 
modelling potentials, e.g., 

V dil = ^{D 2 + D- 2 - 2) = J (ch2g - 1) , k > 0. 
o o 

Obviously, this potential is positively infinite at q = ^foo (D = 0, D = +00) 
and has the global stable equilibrium at q = (D = 1), where it behaves as the 
harmonic oscillator: Vdii(g) « Kq 2 /2 for q « 0. For strongly extended bodies it 
also behaves harmonically in the _D-variable sense. Another phenomenological 
model would be just the global form nq 2 /2. One can also try to use some 
toy models predicting " dissociation" of the body (its unlimited size-expansion) , 
unlimited collapse, or both of them above some threshold of the total dilatational 
energy, e.g., 

VM = ^(\h 2 q-l). 

In certain problems it may be reasonable to use phenomenological models pre- 
venting contraction but admitting dissociation. 

In quantized version of the theory one can stabilize dilatations in an easy 
way with the use of the q- variable potential well (perhaps with the infinite walls) 
concentrated around q = (D = 1). 

If we identify analytically U and V with R n and LI(U,V) with GL(n, M), 
then it is clear that the connected component of unity GL + (n,R) becomes the 
direct product GL+(n,R) ~ SL(n,R) x exp(M) = SL(ra,R) x R+; the second 
group factor is obviously meant in the multiplicative sense, as GL + (1,R). It de- 
scribes pure dilatations, whereas SL(n,R) refers to the isochoric motion. With- 
out this identification, LI([/, V) may be represented as the Cartesian product of 
any of the aforementioned leaves (of mutually non-compressed configurations) 
and the multiplicative group R\{0}. If some volume-standards Uu, vy (e.g., 
metric-based ones v v , v g ) and orientations are fixed in U, V, then LI + (C7, V), 
i.e., the manifold of orientation-preserving isomorphisms, is identified with the 
product SLI([/, vu] V, vy) x exp(R), where, obviously, the first term consists 
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of transformations ip for which A(tp) — 1, i.e., q((p) = 0. Such a formulation 
is more correct from the point of view of geometrical purity, however, for our 
purposes (qualitative discussion of the general solution), the analytical identi- 
fication of U+(U,V) with GL+(n, R) ~ SL(n,R) x exp(M) is sufficient and, as 
a matter of fact, more convenient. In any case, qualitative analysis of the gen- 
eral solution (bounded and non-bounded trajectories) is not then obscured by 
cosmetical aspects of geometry. Thus, from now on the internal configuration 
space Q int = LI(U,V) will be identified with GL+(n, R) ~ SL(n,R) x exp(R). 
Any matrix tp 6 GL + (n,R) will be uniquely represented as (f = Z\P = exp(q)^, 
where ^ € SL(n, R). It is convenient to introduce the following shear velocities: 

d * T -i - T -i<** 
uj := — * \ uj :— * 1 — . 

dt dt 

Obviously u,Cj € SL(n,R)', i.e., they are trace-less. Then affine velocities may 
be expressed as follows: 

= c + ^/, h = u+ d -ll, 
dt dt 

where, obviously, I denotes the identity matrix. 
Analogously, the affine spin splits as follows: 

E = cr + -J, £ = & + -!, cr,a£ SL(n, RV, 
n n 

where p denotes the dilatational canonical momentum. The velocity-momentum 
pairing becomes Tr(EO) = Tr(£f2) = Tr(cra;) + pq — Tr(ercj) +pq. 

Poisson-bracket relations for cr-components are based on the structure con- 
stants of SL(n,R). The same is obviously true for a with the proviso that the 
signs are reversed. The mixed {<r, a} brackets do vanish. Obviously, {q,p} = 1, 
and the quantities q, p (dilatation) Poisson-commute with a, a (shear). 

The doubly-invariant " kinetic energy" PD)) is a superposition of the isochoric 
and dilatational terms, 

A , 9x n(A + nB) , m 
T = -TV) + >-q 2 = T sh + T dil . 

Performing the Legendre transformation, a = Auj, p = n{A + nB)q, we obtain 
the following geodetic Hamiltonian: 

T ^ Tr ^ + M^) P ^ rsh + Tdil - (60) 

In these expressions the quantities u>, a may be replaced by their co- moving 
representations u), a. Lagrangians and Hamiltonians of systems with stabilized 
(controlled) dilatations have the following forms: 

L = L sh + L di] = T sh + T di i - V(q), H = H sh + H dil = T sh + T di] + V(q). 
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There is a complete separability of shear and dilatation degrees of freedom; they 
are mutually independent. This property would not be violated if we included 
also a shear potential V a h,(^), i.e., if V(^,q) = V s h{^) + Vdii(q). The question 
arises as to the structure of general solution for the geodetic SL(n, R)-model, 
i.e., for Vsb — 0. A superficial reasoning based on the analogy with d'Alembert 
models might suggest that the general solution consists only of unbounded mo- 
tions (and the rest states), because there is no potential and the configuration 
space is non-compact. However, it is not the case; there is an open subset 
consisting of bounded orbits. 

Indeed, let us assume that some trace-less matrix a G SL(n, R)' is similar to 
an anti-symmetric matrix A G SO(n, R)', i.e., there exists such \ S SL(n, R) that 
a = x^X 1 - Then every motion ^(t) = cxp(at) 1 i' n is bounded. Indeed, exp(Xt) 
is a bounded subgroup of SO(n,R) C SL(n,R) and so is exp (at) — xexp(Ai)x _1 - 
But similarities are globally defined continuous mappings, therefore, they trans- 
form bounded subsets onto bounded ones. Let us observe that for physical di- 
mensions n = 2, 3 motions of this type are periodic. For higher dimensions peri- 
odicity is not necessary, although obviously possible. To see this, let us consider 
the simplest situation n = 4 and represent R as the direct sum of two comple- 
mentary R 2 -subspaces. Now, let A G SO(n,R)' be a block matrix consisting of 



two skew-symmetric blocks vi 



and V2 



-1 

1 



Rotations gen- 



-1 

1 

erated by separate blocks are obviously periodic, but the total motion is periodic 
if and only if v\jv2 G Q, i.e., the ratio of angular velocities is a rational number. 
If it is irrational, the subgroup obtained by exponentiation of ^-multiples of the 
above block matrix is not a periodic function of the parameter t. It is not a 
closed subset at all; its closure is a two-dimensional submanifold of SO(4,R) C 
SL(4,R). Therefore, being an algebraic subgroup of SO(4,R) C SL(4,R) it is 
not its Lie subgroup in the literal sense. The same concerns any subgroup of 
SL(4,R) obtained from the above one by a similarity transformation. For an 
arbitrary n, solutions of this kind are matrix-valued almost periodic functions 
of the time variable t. 

If a G SL(n,R)' is similar to a symmetric matrix n G SO(n,R)', a = x K X ; 
then, obviously, the motion given by ^(t) = exp(ai)5'o is unbounded. But one 
can show that the previously described bounded almost periodic solutions are 
"stable" in such a sense that for any skew-symmetric A there exists some open 
range of symmetric n-s such that for a = A + k, or, more generally, for similar 
matrices a = x(A + the corresponding solutions ^(t) = exp(crf)'I , o are 

also bounded, although not necessarily almost periodic [43l [74] (and not nec- 
essarily periodic in dimensions n = 2,3). The arbitrariness of pairs (a, ^o) is 
sufficient for the corresponding family of bounded solutions to be open in the 
general solution manifold, thus, 2(n 2 — l)-dimcnsional (topological and differen- 
tial concepts like openness and dimension are meant in the sense of the manifold 
of initial conditions). Let us observe that this statement would be false for solu- 
tions with generators a similar to skew-symmetric matrices. At first look, this 
might seem strange, because the structure of SL(n,R)' implies that the family 
of such as is (n 2 — l)-dimensional and so is the set of initial configurations 
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^o- But these data are not independent and mutually interfere in the formula 
*f?(t) = exp(at)^o- Therefore, the very interesting subfamily of almost-periodic 
solutions is a proper subset of the manifold of all bounded solutions. 

By "anti-analogy", for symmetric matrices A = A T S SL(n,R)' the corre- 
sponding solutions ty(t) = exp(Ai)\&o are non-bounded (escaping in SL(n, R)) 
and, obviously, so are the solutions generated by matrices similar to symmet- 
ric ones, ty(t) — xexp(At)x _1 *o = exp(xAx _1 i)^'o- And again this property 
is stable with respect to small perturbations of A by skew-symmetric matrices 
e = — e T e SO(n,R)'. Therefore, the general solution of the geodetic doubly- 
invariant model contains also an open subset of non-bounded (escaping) solu- 
tions. 

Roughly speaking, using analogy with the Kepler or attractive Coulomb 
problem we may tell here about motions below and above dissociation threshold, 
however, without any potential, just in purely geodetic models on the non- 
compact manifold SL(n, R). 

In a quantized version of this model the family of bounded classical trajecto- 
ries is replaced by the discrete energy spectrum and the L 2 -class wave functions 
of stationary states. And similarly, the manifold of non-bounded orbits is a clas- 
sical counterpart of the continuous spectrum and non-normalized wave functions 
(scattering situations). 

Obviously, the above description in terms of groups GL(n,R), SL(n, R) 
is an analytical simplification used for computational purposes. To use sys- 
tematically a more correct geometrical language we should replace the terms 
"skew-symmetric" and "symmetric" by g- or ^-skew-symmetric and symmetric: 

\*j = T9 lk X m kgrnj, X A B = TV AC ^ D CVDB- 

Finally, we can conclude that if dilatations are stabilized by some poten- 
tial Vdn(q), then for the model with the kinetic energy invariant under spatial 
and material affine transformations, there exists a 2n 2 -dimensional family of 
bounded solutions even if the shear component of motion is purely geodetic. 
If the stabilizing dilatational potential has an upper bound, there exists also 
a 2n 2 -dimensional family of unbounded, escaping motions. The above argu- 
ments are based on properties of one-parameter subgroups and their cosets in 
SL(n, R). Therefore, they do not apply directly to afFine-metrical and metrical- 
affine models. Indeed, as we have seen, if the spatial or material symmetry of 
the kinetic energy is restricted to the rotation group, then, except some special 
solutions (relative equilibria), one-parameter subgroups and their cosets fail to 
be solutions. Nevertheless, our arguments may be used then in a non-direct 
way. 

In analogy to (|60|) we can rewrite the kinetic Hamiltonians (|4"Tj) . ([42]) . i.e., 
[(IS ]) , (|3S] 1 . as follows: 

^ = WTAj^ + 2n(I + A + n B / + 2(7^) (61) 

1 -c SHn) (2)+ —- T . \ . , p y + ' A2 J V\\ 2 , 



2(1 + A) ou{n,y ' 2n(I + A + nBy 2(I 2 - A 2 ) 
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2(1 + A) °^ n,y ' 2n(I + A + nBy 2(I 2 - A 2 ) 

where C SL(n) (2) = Tr(a 2 ) = Tr(a 2 ). 

The formulas (|45|) , (|46|) or, equivalently, (|6Tj) , ([62]) imply that for dilatatio- 
nally-stabilizcd models H = 7I n t + Vdn{q) with the affine- metrical and metrical- 
afhnc kinetic terms % n t, all the above statements concerning bounded and un- 
bounded solutions of afhne-afhne models (fBTTj) . (|2T))) remain true. In particular, 
for the purely geodetic incompressible models with % nt invariant under SL(V^)x 
0(U,rj) or under 0(V,g)x SL({7), there exists an open subset of bounded so- 
lutions (vibrations) and an open subset of non-bounded ones. What concerns 
spatially affine and materially metrical models, the very rough argument is that 
the evolution of quantities S, 1C is exactly the same as it was for Hamiltonians 
H with 7i nt affinely-invariant both in the physical and in the material spaces, 
in this case in (|60[) A is replaced by I + A. This is a direct consequence of 
equations of motion written in terms of Poisson brackets, 

In fact, \\V\\ is a constant of motion for both types of Hamiltonians (affine-affine 
and affine- metrical) . In addition to Lie-algebraic relations of GL(V)' ~ LfV) 
satisfied by we have the following obvious Poisson rules: 

C(2)} = {£*,-, C(1)} = 0, {E i j ,\\V\\ 2 } = 0, {!C a ,\\V\\ 2 } = 0. 

The first equations express an obvious property of C(k) as Casimir invariants 
of T, l j (and Y> A b)- The second formula follows from the obvious relationship 
{E'jjS^s} = 0, because ||V|| 2 is an algebraic function of Y* A b- And the third 
equation is due to the fact that the deformation invariants IC a are invariant 
under the group of material isometries generated by V A b- 

Therefore, the time evolution of variables JC a is identical in both types 
of models, i.e., (|60p and (16 1 [) ; the former with A replaced by / + A. As a 
matter of fact, for geodetic models with dilatation-stabilizing potentials V(q), 
the deviator a l j = S l j — (1 / 'n)S a Q <5 1 j is a constant of motion and, obviously, it 
is so for the purely geodetic incompressible models. The only difference occurs 
in degrees of freedom ruled by SO(V,g), SO(U,r}), describing the orientation of 
principal axes of deformation tensors C € V* <E> V*, G £ U* <£>U*. But, roughly 
speaking, these degrees of freedom have compact topology and their evolution 
does not influence the bounded or non-bounded character of the total orbits. 

The same reasoning applies to dilatationally stabilized geodetic models in- 
variant under 0(V,g) X GL(U) or purely geodetic incompressible models invari- 
ant under 0(V,g) x SL((7) (spatially metrical and materially affine models). 
Then everything follows from Poisson brackets 

{£ A B ,C(2)} = {£ A B ,C(1)} = 0, {± A B ,\\S\\ 2 } = 0, {/C a ,||5|| 2 } = 0. 
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Now on the level of state variables £ A b, K, a the time evolution is exactly iden- 
tical with that based on the afnne-afnne model of % nt (again with A in (|60|) 
replaced by I + A). 

Let us stress an important point that it is the time evolution of deformation 
invariants that decides whether the total motion is bounded or not. This is a 
purely geometric fact independent on any particular dynamical model. There 
is an analogy with the material point motion in K™. An orbit is bounded if and 
only if the range of the radial variable r is bounded. 

The above point plays an essential role in the qualitative discussion of de- 
formative motion. It suggests one to use analytical descriptions of degrees of 
freedom based on deformation invariants. 

6 Analytical description 

In section 3 some fundamental facts concerning deformation tensors and defor- 
mation invariants were summarized. Below we continue this subject and present 
some natural descriptions of affine degrees of freedom well-adapted to the study 
of isotropic problems. 

The material and physical spaces are endowed with fixed metric tensors, 
ij £ U* ®U* , g € V* ® V*, and any configuration ip € LI(£/, V) gives rise to the 
symmetric positively definite tensors G[<p] £ U*®U*, C\ip] £ V*<S>V*, i.e., Green 
and Cauchy deformation tensors. Raising their first indices respectively with 
the help of rj and g, we obtain the mixed tensors G[ip] £ U (£> U* , C[tp] £ V <g> V* 
with eigenvalues A„, A~ , o = 1, n. It is also convenient to use the quantities Q a , 
q a , where Q a = exp(q a ) = y/X^. The diagonal matrix D — diag(Q 1 , . . . ,Q n ) is 
identified with the linear mapping D : M" — > R n . 

The configuration <p £ LI(J7, V) may be characterized by D, i.e., by the 
system of fundamental stretchings Q a — exp(q a ), and by the systems of eigen- 
vectors R a £ U, L a £ V of G, C normalized, respectively, in the sense of r\ and 
9, 

GR a = X a Ra = exp(2 g a )i? Q , CL a = \- x L a = exp(-2<f )£„. 

Obviously, when the spectrum is non-degenerate, then R a , L a are uniquely de- 
fined (up to re-ordering) and pair-wise orthogonal, r](R a , Rb) = t/cdR aR D b — 
$ab = g(L ai L a ) = gijU a L\. Such a situation is generic, thus, when at some 
time instant (Gl (fi(t) corresponds to degenerate situation, then L a (t), R a (t) 
may be also uniquely defined due to the continuity demand. 

The elements of the corresponding dual bases will be denoted respectively 
by R a £ U* , L a £ V* . When necessary, to avoid misunderstandings, we shall 
indicate explicitly the dependence of the above quantities on cp £ L1(U, V): 
q a [ip], R a [ip], L a [ip], etc. 

Green and Cauchy deformation tensors may be respectively expressed as 
follows: 

G[<p] =X>«M#>]®#>] =X>p(2g°M)ii>]®i? a M, 

a a 
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CM = E A aV™ ® = E 6X P ("VM) £ » ® L< M 

a a 

In this way <p has been identified with the triple of fictitious objects: two 
rigid bodies in U and V with configurations represented, respectively, by or- 
thonormal frames R G F(U, 77), L £ F(V, 5) and a one-dimensional n-particle 
system with coordinates q a (or Q a ). Even for non-degenerate spectra of G[ip], 
C[tp] this representation is not unique because the labels a under the summation 
signs may be simultaneously permuted without affecting ip itself. For degener- 
ate spectra this representation becomes continuously non-unique in a similar 
(although much stronger) way as, e.g., spherical coordinates at r = 0. 

Let us observe that the linear frames L = (. . . , L ai . . .) and R = (. . . , R a , . . .) 
may be, as usual, identified with linear isomorphisms L : R™ — > V and R : R n — > 
Z7. Similarly, their dual co-frames L = (. . . , L a , . . .) and R = (. . . , i? a , . . .) are 
equivalent to isomorphisms L^ 1 : V — > K n and : CI — * R n . Identifying the 
diagonal matrix diag(. . . , Q a , . . .) with a linear isomorphism D : R™ — > R n , we 
may finally represent 

if = LDR' 1 , 

this is a geometric description of what is sometimes referred to as the two-polar 
decomposition [H EH E3 [Ml E3 EE] . 

Strictly speaking, in continuum mechanics, when the orientation of the body 
is constant during any admissible motion (no mirror- reflections), one has to 
fix some pattern orientations in U, V and admit only orientation-preserving 
mappings tp. And then the non-connected sets of all orthonormal frames F(£/, rj), 
F(V, g) are to be replaced by their connected submanifolds F + ([7, 77), F + (V,g) 
of positively oriented frames. 

Obviously, the spatial and material orientation-preserving isometries A £ 
SO(V,g), B S SO(U, 77) affect only the L- and i?-gyroscopes on the left. Indeed, 
L 1— > AL, R 1— > BR result in ip 1— > A(^_B _1 . Their Hamiltonian generators, spin 
and minus- vorticity (i.e., respectively V- and [/-spatial canonical spins) have 
identical Poisson-commutation rules. 

For any of the mentioned rigid bodies, one can define in the usual way the 
angular velocity in two representations. One should stress that both V and U are 
from this point of view interpreted as "physical spaces". The "material" ones 
are both identified with K™. The "co-moving" and "current" representations 
X G SO(n, R)', x € SO(V,g)' for the L-top are respectively given by 

X 5 :=< L , — >= L t —, X := X b L a ® L , ue., ^ = 

The corresponding objects 7) G SO(n, R)', 7? G SO(U, rj) 1 for the i?-top are 
defined by analogous formulas: 

d\ :=< R a , ^ >= R a K^, := # a b R a ®R b , i.e., - ^^l- 

In certain problems it is convenient to use non-holonomic velocities g a , x a fc, ?9 a b 
or q a , x % j-, "& A b- Similarly, non-holonomic conjugate momenta p a , p a b, T a b or p a , 
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p l j, t a b are used, where again p,f e SO(n,R)', p € SO(V,#)', r e SO(f7,7?)'. 
The pairing between non-holonomic momenta and velocities is given by 

< (p,T,p),(x,#,q) >=< {p,T,p),{x,fi,q) > 

= Pa q a + ^Tr(px) + ^Tr(rtf) = p a q a + ^Tr(px) + ^Tr(ftf). 

Remark: Our system of notations is slightly redundant, because p and r coin- 
cide, respectively, with spin and negative vorticity, p = S, r = —V. The reason 
is that they arc Hamiltonian generators of transformations ip i— » Aip, ip i— > ipB^ 1 , 
AeSO(V,g), BeSO{U,r)). 

The objects p, f generate transformations 

L^LA, R^RB, A,BeSO(n,R) (63) 

and express the quantities p, r in terms of the reference frames given, respec- 
tively, by the principal axes of the Cauchy and Green deformation tensors, 

P = p a b L a ®L b , t = r a b R a O R b . 

Remark: In dynamical models based on the d'Alembert principle the quantities 
Q a and their conjugate momenta P a are more convenient that q a and p a . The 
latter ones are useful in models with affinely-invariant kinetic energy. 

If V and U both are identified with R n and LI(U, V) with GL(n, M), then L 
and R in the two-polar splitting <p = LDR' 1 become elements of SO(n, R) and 
D, as previously, is a diagonal matrix with positive elements. The two-polar 
decomposition is a by-product of the polar decomposition of GL+(n,R), 

<p = UA, 

where U € SO(n,R), thus, U T — U^ 1 1 and A — A T is a symmetric positively- 
definite matrix. It is well-known that this decomposition is unique, whereas 
the two-polar one is charged with some multivaluedncss. Green and Cauchy 
deformation tensors are then represented as follows: G = (fi T <p = A 2 , C = 
(ip~ 1 ) T ip~ 1 = ?7A _2 f _1 . One can also use the reversed polar decomposition 

(p = BU, Ue SO(n, R), B = UAU' 1 = B T . 

Then G = U~ 1 B 2 U, C = B~ 2 . The two-polar decomposition is achieved by the 
orthogonal diagonalization of the matrix A, A = VDV^ 1 , V € SO(n,R). Then 
L = UV, R = V. 

The polar splitting was described above in an over-simplified standard way, 
namely, U and V were identified with R" and LI(U,V) with GL(n, R). Let 
us remind that in continuum mechanics the connected components of LI([/, V) 
and GL(n, R) are used as configuration spaces, LI + ([7, V), GL + (n,R), where 
the first symbol denotes the manifold of orientation-preserving isomorphisms 
(it is assumed here that some orientations in U, V are fixed). It is instructive to 
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see what the both polar splittings are from the geometric point of view, when 
U and V are distinct linear spaces, non-identified with R™. 

As mentioned above, when metric tensors rj £ U* <g> U* , g £ V* ® are 
fixed, then any <p £ LI(C7, V) with non-degenerate spectra of deformation tensors 
gives rise to the pair of orthonormal bases (L a [ip] £ V, a = l,n), (Ralf] £ 
U, a — l,n). There exists exactly one isometry {/[y] : U — > V such that 
U[ip] ■ Raltp] = LoM- Obviously, the isometry property is meant in the sense 
that r) = U[(fi]*-g, i.e., analytically tjab — gijUYp] 1 aU[lp\^ b ■ Geometric meaning 
of the polar decomposition is as follows: 

<p = U[ip]A[(p] = B[tp]U[<p], 

where the automorphisms A[(p] £ GL(U), B[ip] £ GL(V) are symmetric, re- 
spectively, in the rj- and g-sense, i.e., r)(A[(p]x,y) = r)(x, A[tp]y), g(B[ip]w, z) = 
g(w, B[tp]z) for arbitrary x,y £ U , w, z £ V . They are also positively definite, 
r](A[tp]x, x) > 0, g(B[tp]w, w) > for arbitrary non-null x £ U, w £V. 

In spite of the non-uniqueness contained in L[ip], R[f], the mappings U[tp], 
A[ip], B[(p] are unique. And the symmetric parts are obtained from each other 
by the J7[<^]-intertwining, B[ip] — f7[<^]j4[<^]J7[<^] _1 . 

In mechanics of discrete affinc systems we are free to admit orientation- 
reversing isometries U or symmetric mappings A, B not necessarily positively- 
definite. 

The non-uniqueness of the two-polar decomposition mentioned above is im- 
portant in certain computational and also principal problems, so some comments 
are necessary here. The problem is technically complicated, thus, only necessary 
facts will be quoted here, some of them formulated in a rather brief, rough way. 

The subgroup of 0(n, R) consisting of matrices which have exactly one non- 
vanishing entry in every row and column will be denoted by K. Obviously, K 
is finite and the mentioned entries are ±1, reals with absolute value 1. The 
subgroup of proper if -rotations will denoted by K + := K (~l SO(n,R). Ob- 
viously, the orders (numbers of elements) of K, K + equal respectively 2n ■ n\ 
and n ■ n\. Let W £ K be a corresponding similarity transformation preserving 
the group of diagonal matrices, Diag(R") 3 D W~ 1 DW £ Diag(R"), and 
resulting in permutation of the diagonal elements of D — diag(Q 1 , . . . , Q n ), i.e., 
we have (Q\...,Q™) ^ (Q 7Tw{1 , . . .,Q n) ) or (q 1 ,...,q n ) h-> (g^ 1 , . . . , q n )), 
where Q a = exp(q a ). Obviously, the mapping K 3 W ^ ttw £ is a 2n : 1 
cpimorphism of K onto the permutation group S( n ). Its restriction to K + has 
an n-element kernel. The non- uniqueness of representation of ip £ GL + (n, R) 
through elements of SO(n,R) xl"x SO(n, R) depends strongly on the degen- 
eracy of spectra of deformation tensors. The multi-valuedness is discrete, thus, 
simplest in the case of simple spectra. 

Let GL+(™)(ra,R) C GL + (n,R) be the subset of ip-s with non-degenerate 
spectra of C, G. The corresponding subset AfW of SO(n,R) x R" x SO(n,R) 
consists of such triplets (L; q 1 , . . . ,q n ; R) that all q l -s arc pairwise distinct. The 
group K + may be faithfully realized by the following transformation group 
ofMW: 

(L; q 1 ,..., q n ; R) ~ (LW; q^\ . . . , q^-RW). 
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Obviously, this transformation does not affect tp = LDR , Therefore, we have 
a diffconiorphism 

GL +(n) (n,IR) ~M(")/ff( n ). 

Non-degenerate spectrum is a generic one, nevertheless the coincidence case 
must be also taken into account because some new qualities appear then and 
they are relevant for qualitative analysis of classical phase portraits and for 
quantum conditions on admissible wave functions. 

Let GL + ( k ' Pl >'"' Pk ' C GL + (n,R) consist of ip-s for which deformation tensors 
have k < n different principal values, every one of them with the corresponding 
multiplicity p a , Y^ k a=1 Pa = n. And similarly, let M( k ' ,pl '"' ,Pk> be the set of 
such triplets (L; q 1 , . . . , q n \ R) S SO(n, R)xM"x SO(n, R) that there are only k 
different q l -s with the same conditions concerning multiplicity. And now let us 
consider the transformation group H^ k ' pi ^—' Pk > acting on M^ k ' pl '"' ,Pk " > as follows: 

(L; q\ . . . , q n -R) -> (LW; . . . ; q n )-RW), 

where W runs over the subgroup of SO(n, R) that is generated by K + and the 
subgroup 

H (k; Pl ,..., Pk ) c sO(n,R) composed of k blocks Pa x Per, every one given 
by the corresponding SO(p ff) M). Then we have that 

Qj_+(k;pi, ■■-,Pk) ^ M(k;pi,-~,Pk) I jj{k;px,---,Pk) _ 

When k < n, then at least one of multiplicities is non-trivial and the resulting 
group H^ k '' Pl ''"' Pk ' is continuous. The resulting quotient is lower-dimensional 
because of this continuity of the divisor transformation group. 

In the physical case n = 3, we have obviously only two possibilities of the 
non-trivial blocks, namely the total SO(3,R) and SO(2,R) x SO(l,R) (respec- 
tively, all three q l, s equal or two of them); obviously SO(l,R) = {1}. 

In the extreme case k = 1, D is proportional to the n x n identity matrix and 
it is only the total LR~ X that is well-defined; on the other hand L, R separately 
are meaningless. 

It is very convenient and instructive to express our Hamiltonians, kinetic 
energies and configuration metrics in terms of the two-polar splitting. The 
previous statements concerning the phase pictures become then much more 
lucid. Let us introduce some auxiliary quantities M := — p — f , N := p — f. 
One can easily show that the second-order Casimir invariant C(2) occurring in 
the main terms of our affine-affine, affine-metrical and metrical-affine kinetic 
Hamiltonians has the following form: 

c ( 2 ) = 2^ + i^^z^-i^5?Sr (64) 

a a.b hlL 2 a.b Ln 2 

Obviously, M and N are antisymmetric in the Kronecker-delta sense, M a b = 
-M b a = -g bk g al M k h N a b = -N b a = -g bk g al N k i. The first term in flM|) may 
be suggestively decomposed into the "relative" and the "over-all" ("centre of 
mass") parts: 

1 . 

Zn * — ' n 

a,b 
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Obviously, C(l) = p. 

For geodetic systems and for more general systems with potentials V de- 
pending only on deformation invariants, spin S — p and vorticity V = r are 
constants of motion and may be used for extracting from equations of motion 
some information concerning the general solution. Unlike this the quantities 
p, f , thus, also M, N, fail to be constants of motion except the special case 
n — 2, when the rotation group is Abelian. However, on the level of qualita- 
tive analysis, the expression ([SI]) based on p, f is more convenient because it 
does not involve L, i?-variables, i.e., rotational degrees of freedom of deforma- 
tion tensors. Therefore, our Poisson bracket relations imply that on the level 
of variables q a , p a , M a b, N a b equations of motion based on ([ST]) (equivalently 
©), © (equivalently gBJ), and (JBDJ) with A replaced by I + A are identical. 
In particular, for geodetic incompressible models and for compressible models 
with stabilized dilatations there exists an open family of bounded (vibrating) 
solutions and an open family of non-bounded (decaying) solutions. The reason 
is that it is so for ([SO)) with A replaced by I + A, and the additional terms 
proportional to S 2 or V 2 do not influence anything because they have vanishing 
Poisson brackets with q a , p a , M a i,, N a i, and only those variables occur in H. 
The only difference appears when the evolution of L- and i?-variables is taken 
into account. However, the corresponding configuration spaces F(V,g), F(U,r)) 
are compact (they are manifolds of orthonormal frames) and do not influence 
the boundedness of orbits. 

Let us observe that after substituting ([64]) . the first main term of ([ST]) (equiv- 
alently (J45])), ((62]) (equivalently ([IS])), and (J60J with A replaced by I+A acquires 
the characteristic lattice structure, 

= J_ v 2 1 v (M\) 2 1 V (N\) 2 

latt 2a 2^ Pa + 2^ ,2 q °-a» 32 a ^ ch 2 <£=<t ' 

a a ,b ° 1L 2 a,b L11 2 

This expression resembles structurally the hyperbolic Sutherland n-body sys- 
tem on the straight line. Positions of the fictitious material points are given by 
deformation invariants q a . The "particles" have identical masses and are indis- 
tinguishable. Unlike the hyperbolic Sutherland system, the coupling amplitudes 
M a b, N a (, are non-equal and non-constant; rather they are dynamical variables 
on the equal footing with q a , p a . The negative A^-contribution to T\ &tt describes 
the attractive forces between lattice points, whereas the positive M-term cor- 
responds to repulsion. Under the appropriate initial conditions we have stable 
bounded vibrations without any use of the potential energy term. Therefore, 
the non-definiteness of 7i a tt is not only non-embarrassing, but just desirable as 
a tool for describing "elastic" vibrations on the basis of purely geodetic models. 
Let us observe that the purely affine-affine part of ([ST]) , ([S2"]) (equivalently (|4"5]) , 
(J3SJ), i.e., ([SO]) with A replaced by I + A (composed of its first two Casimir 
terms), splits in the following suggestive way into the binary SL(n,R)-part and 
dilatational contribution: 

a,b 
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(M a b ) 2 1 ^ {N\) 2 , na + /3 2 

P > 



32a^ sh 2 £^ 32a ^ C h 2 £^ 2na/3 



or, in a more explicit form, 

/-raff 1 \~"V \2 1 {M a 



32(7 + A) ^ ch 2 2l_V 2n(I + A + nB) P ' 1 j 

Obviously, for (fBTj) (equivalently (|45p) and (f6"2")l (equivalently (|46p ) we have, 
respectively, 

^f- m ° tr = + 2{p [ A2) \\Vf, (66) 

^etr-aff = ^aff + _ L_ _ | , g| , 2 (g7) 

Copmparing this with ([6"T]) . (|62p . we conclude that 

C SL(rl) (2)=Tr(a 2 ) = Tr(a 2 ) 

- 1 \^ _ ^ , 1 V- ( MO fc) 2 1 v-JA[V_ 

- 2n ^ a Pbj + 16 sh 2£^L 16 2- ch 2 «^ ' 

a,fc a,b B11 2 a,6 2 

This expression is very suggestive because it expresses the quantity C SL („)(2) 
and the corresponding contribution to 7i nt , i.e., the metric tensor on the mani- 
fold of incompressible motions, as the sum of n(n ~ l)/2 two-dimensional clus- 
ters, i.e., R 2 -coordinate planes in R n . Incompressibility is expressed by the fact 
that the invariants q a and their conjugate momenta p a enter the above formula 
through the shape-describing differences (q a — q b ) (ratios Q a /Q b ) and p a — Pb- 
This expression may be very convenient when studying invariant geodetic mod- 
els on the projective group Pr(n,R), i.e., when dealing with the mechanics of 
projectively-rigid bodies (bodies subject to such constraints that all geometric 
relationships of projective geometry arc preserved, in particular, the material 
straight lines remain straight lines). The point is that Pr(n, R) may be identified 
in a standard way with SL(n + 1, R). 

For the d'Alembert model the two-polar splitting leads to the following ki- 
netic Hamiltonian term: 

T. 1 V P 2 1 V ( Ma ") 2 , 1 V ( Na rf (6R ) 
mt 21 ^ a 81 {Q a - Q b ) 2 81 (Q a + Q b ) 2 ' 1 ' 

It is purely repulsive on the level of Q-variables, thus, without any potential 
term it is non-realistic as a model of elastic vibrations. It is related to the 
Calogero-Moser lattices similarly as the previous models show some kinship 
with the hyperbolic Sutherland lattices [iSDHHHSllSllMllMlED]. 
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What concerns affinc models, we can compactify deformation invariants q a 
by taking them modulo 2ir (n-dimensional torus), i.e., by putting formally 
Q a = cxp(iq a ). This is equivalent to replacing GL(n, R) by U(n), i.e., an- 
other and completely opposite real form of GL(n, C). The Lie algebra U(n)' 
consists of anti-Hermitian matrices, and the positively definite kinetic energy 
may be postulated in the following form: 

T in t - -^Tr(r> 2 ) - |(Trf!) 2 = ^Tr(f!+0) + |Tr(Q+)Tr(n), 

where = (dip/dt)^ 1 , A > 0, B > 0. Obviously, in this expression for T, £1 
may be as well replaced by f2 = ip -1 (dip / dt) . 

Using again the "two-polar" decomposition ip = LDR^ 1 , where L,R e 
SO(n, R), D = diag(. . . , exp(ig a ), . . .), one obtains for the geodetic Hamiltonian: 

r t - 1 W i 1 V (Ma&)2 i 1 V {Na » )2 * p 2 

mt 2A^ Pa i2A^ s ^^t 2,2A ^ cos 2 sl^t 2A(A + nB) 

o, (2,6 2 ci,6 2 

(69) 

The first three terms, corresponding to the C(2)-Casimir, resemble the usual 
Sutherland lattice for g-particles with the same provisos as previously. Geodetic 
motion is bounded, because U(n) is compact. Just as previously, it may be 
convenient to use the splitting into SU(n)- and U(l)-terms, 

^ = TAa E(P- - + 2n(A + n B / 

a.b 



32^ ^ sin 2 32A ^ 

a, 6 2 a,b 

And, in particular 



(iV a 



\2 



C su(n )(2)--^ a -p 6 ) + _ 



a,b - a,6 Sm 2~ ±U a.fc COS^ 



The binary structures of C,g£(n,R)(2) and Csj/(„)(2) and their dependence 
on the variables q a , p a through their differences q a — q b , p a — pb is geometrically 
interesting in itself. The splitting into SL(2, R)- and SU (2)-clusters correspond- 
ing to all possible coordinate planes R 2 in R" may be also analytically helpful. 
However, some sophisticated mathematical techniques would be necessary then, 
like, e.g., the Dirac procedure for degenerate/constrained system. The point 
is that, in general, different clusters arc not analytically independent. And 
any procedure based on some ordering of variables destroys the explicit binary 
structure and makes the structure of T rather obscure. 

It is interesting that the general solution of C(2)-based geodetic models con- 
tains as a particular subfamily the general solution of the mentioned Calogero- 
Moser and Sutherland models. It is obtained by putting N a b = 0, and all M a b 
with b ^ a equal to some fixed constant M. 
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As mentioned, we are particularly interested in geodetic affine models. Nev- 
ertheless, it is instructive to admit a wider class of Hamiltonians: 



H = T + V(q\...,q n ), (70) 

where T is any of the kinetic energy models described above, and the potential 
V depends on (p only through the deformation invariants q a . This means that 
it is isotropic both in the physical and material space. The mentioned non- 
uniqueness of the two-polar decomposition implies that V as a function on 
R™ must be permutation-invariant to represent a well-defined function on the 
configuration space. When the extra potential, e.g., elastic one, is admitted, 
then also the "usual" model of T based on the d'Alembert principle may be 
sensibly used for describing bounded elastic vibrations. Therefore, from now 
on all the above models of T int (fu5)l . (fuT)) . are admitted, although 

of course the "non- usual" affine models (p)5|) , (ffiu]) . (pjTj) are still particularly 
interesting for us. 

As mentioned, the most convenient way of discussing and solving equations 
of motion is that based on Poisson brackets, 

where F runs over some maximal system of (functionally) independent functions 
on the phase space. The most convenient and geometrically distinguished choice 
is q a , p a , M a b, N a , L, R or, more precisely, some coordinates on SO(n,R) 
parameterizing L and R. In d'Alembert models Q a , P a are more convenient the 
Q a , Pa- 

An important point is that q a , p a , M a , N a generate some Poisson subal- 
gebra, because their Poisson brackets may be expressed by them alone without 
any use of L, R- variables. And Hamiltonians also depend only on q a , p a , M a b, 
N a , whereas L, R are non-holonomically cyclic variables. This enables one to 
perform a partial reduction of the problem. In fact, the following subsystem of 
equations is closed: 

dqa i# m dH 

at dp a 
= {M\, H} = {M\, M c d }^— + {M%, N c d }- ° 



dt J UJ 9M c d " J cW c d ' 

dpa r dH 

— = {Pa,H} = - — , 

Obviously, {q a lPb } = 6\, {q a ,M c d } = { Pa ,M c d } = {q a ,N c d } = { Pa ,N c d } = 
0. Poisson brackets of M, iV-quantities follow directly from those for p, f , and 
the latter ones correspond exactly to the structure constants of SO(n, R), thus, 

{pab, Pcd} = PadScb — PcbSad + pdbSac — PacSclb, 
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{fab, f cd } — fadScb ™ fcb&ad + fdb$ac ~ fac^db, {Pab,f C d} = 0, 

where the raising and lowering of indices are meant in the Kronecker-delta sense. 
From these Poisson brackets we obtain the following ones: 

{Mob, Mod} = {Nab, N cd } = McbSad - MadScb + M d bS a c ~ M ac S d b, 
{M a b, N cd } = N cb S a d - NadScb ~ N ac S d b + N do 8 a c- 

The subsystem for {q a ,p a , M a b, N a b) may be in principle autonomously solv- 
able. When the time dependence of p = (N - M)/2 and f = -(N + M)/2 is 
known, then performing the inverse Legendre transformation we can obtain the 
time dependence of angular velocities x, 

v» - dH &° dH 

X b - stj-, v b 



dp h a ° df\ 

(some care must be taken when differentiating with respect to skew-symmetric 
matrices). And finally the evolution of L, R is given by the following time- 
dependent systems: 

dL dR * 

There is some very important consequence of this reduction procedure, i.e., 
in doubly-isotropic models spin S, vorticity V, and their magnitudes \\S\\, \\V\\ 
are constants of motion. Moreover, ||5|| and ||V|| have vanishing Poisson brack- 
ets with all quantities q a , p a , M a b, N a b- Therefore, on the level of these vari- 
ables, all Hamiltonian systems (|70p with the same doubly isotropic potential V 
and with three affine models of the kinetic energy ([B5"j). (|57|) are identical. 
In particular, the solutions for variables q a , p a , M a , N a b coincide with those 
for the affine-affine model (|65|) . And this applies, in particular, to the geode- 
tic model (when V = 0) and to the geodetic shear model with extra imposed 
dilatations stabilized by Vdn(q), where q = (q 1 + ■ ■ ■ + q n )/n. And then, as men- 
tioned, the argument about one-parameter subgroups and their cosets decides 
about the existence of open subsets of bounded and non-bounded trajectories. 
The only difference between various T-models appears only on the level of L, 
i?-degrees of freedom. But the compactness of the corresponding configuration 
spaces F(V,g), F(U,rf) implies that this part of motion does not influence the 
property of the total orbits in Q = Ll(U, V) to be bounded or non-bounded. 
One should stress that for the affine-metrical and metrical-affine geodetic mod- 
els ([SB]) . (1571) only exceptional solutions are given by one-parameter subgroups 
and their cosets (relative equilibria). Nevertheless, extracting from all possible 
one-parameter subgroups and their cosets their (q a , p a , M a b, TV^-content, we 
obtain true statements concerning all three geodetic models (|65p. ([55]) . (|67|) . 

Our affine geodetic models (p5|) . ([57)1 have a nice binary structure 

with an additional degree of freedom related to the motion of the centre q 
of deformation invariants q a , a = l,n. In practical applications this term in 
7i nt should be stabilized by some extra introduced dilatational potential. If we 
perturb geodetic models by admitting more general doubly-isotropic potentials, 
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then it follows from the mentioned structure of % nt that the most natural and 
computationally effective potentials will be those somehow adapted to the above 
splitting into shear and dilatation parts, i.e., 

V(q\ . . . , q n ) = V da (q) + ^ V sh ^ (\q* - q>\). 

Here the additional shear part is not only binary but, just as it should be, it 
is depending only on the relative positions of deformation invariants \q l — q 3 \ 
on R or on the circle U(l) when the group U(n) is used. Obviously, the model 
of V^h = (1/2) J2i j Vsh^dl 1 ~~ O'D wm b e computationally effective only when 
the structure of functions will have something to do with sh — q^/2), 
ch ([^-^]/2), sin ([ q i- q i)/2), 008(^-^/2). 
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